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ABSTRACT 



A new method for absolute calibration of electroacoustic trans- 
ducers is described which permits low frequency reciprocity calibration 
in a test apparatus of small dimensions. The method takes advantage of 
the reduced wave speed in a cylindrical water column bounded by a 
compliant PVC tube to reduce the length of the standing wave to 
dimensions which can be easily handled in the laboratory. 

- Results of a limited series of experiments show a mean reproduci- 
bility of iO.8 decibels in a one-meter long tube at frequencies between 
750 and 1100 Hz and sound speeds from 330 to 365 meters/second . The 
hydrophone voltage sensitivities obtained by the method, however, 
systematically differ by 1.5 decibels from those measured by comparison 
with a standard hydrophone; additional research will be required to 
resolve this discrepancy. 
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I. WHV COMPLIANT TUBE RECIPROCITY WORKS 



A. ACOUSTIC CALIBRATION AND RECIPROCITY 

There are two general classes of underwater acoustic transducer 
calibration: absolute, or primary, calibration, and comparison, or 

secondary, calibration. The latter presumes the availability of a sound 
field of known amplitude measured by a calibrated transducer, but 
absolute calibration can be performed using only voltage measurements 
and primary measurements of the test system environment in a sound 
field of unknown amplitudo. Once the voltage sensitivity of a standard 
hydrophone has been established in a particular frequency range and 
environment by absolute calibration, it may be used as a referenc.e for 
comparison calibration of other hydrophones and of projectors; this 
absolute calibration is most often done by the reciprocity method. 

1 . Fundamentals of Reciprocity Calibration 

The reciprocity method for absolute calibration is based on the 
properties of linear, passive systems whose action can be described by 
the two-port equation (1.1), which relates the acoustic pressure and 
velocity (p,u) at a transducer’s mechanical terminal and the voltage and 
current (E,l) at its electrical terminals by a symmetric or antisymmetric 
matrix of constant coefficients; such a system is called reciprocal. 



■ Zn, 

,±<l> Ze J 



( 1 . 1 ) 



Hunt [Ref. 1] derives this equation from first principles and shows how 
the matrix elements may be evaluated for common transducers. Foldy 
and Primakoff show [Ref. 2] that the short circuit transmitting current 
response S and the open circuit hydrophone voltage response M of a 
reciprocal transducer T are related by a constant factor 
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where J is called the reciprocity parameter and is characteristic of the 
test medium and its boundaries. Bobber has demonstrated [Ref. 3] 
that J is equivalent to the acoustic transfer admittance between an 
arbitrarily defined transmitting surface T generating a volume velocity 
U and an arbitrarily defined receiving surface H within the test envi- 
ronment exposed to an average pressure P, and exists as a function of 
the test medium and its boundaries for any arbitrary linear, passive, 
and reversible medium. 

When two test transducers^ are located at points in the 
medium where the pressure amplitudes are equal, the reciprocity param- 
eter allows expression of the signal voltage of a test hydrophone H 
placed at one point as a function of the transmitter current response of 
a reversible transducer T placed at the other point and driven to 
create the sound field; 



Eth — MhPh — MhStIt 



(1.3) 



where I is the drive current input to T. If the sound field generated 
by T is then replaced by a field which is identical at the locations of T 
and H, but is generated by an separate projector P, the outputs of T 
and H in this new field will be related by equation (1.4). Combining 
this equation with (1.3) produces the calibration equation (1.5), which 
allows the computation of the test hydrophone's voltage sensitivity H 
directly from voltage measurements made in the three sou rce-to- receiver 
configurations: T-to-H, P-to-H, and P-to-T. 



^To avoid confusion, the term test transducers will always be used 
when referring to both transducers used in the reciprocity^ calibration . 
The term reversible transducer, or the ANSI designator ’’T", will be 
used to refer to the reciprocal transducer used both as a hydrophone 
and as a projector during the calibration; test hydrophone, or ''H" will 
designate the transducer used only as a hydrophone; and projector, or 
"P" will identify a third sound source used to create the sound field 
used in the test. 
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(1.4) 



EpT - ^tP _ J S, 
Ep„ M„P . M„ 

2 EE 

Mh = J 

EPT It 



(1.5) 



2. Reciprocity Test Geometries 

Most reciprocity calibrations are made in one of the two 
earliest system geometries for which the reciprocity parameter was 
computed: free field [Ref. 4], or close-coupled [Ref. 5]. The tech- 
niques for carrying out these tests, and the limitations required to 
ensure the validity of the test are described by the ANSI standard for 
underwater calibration [Ref. 6]. The free field test geometry requires 
a volume of water large enough to permit spherical spreading of the 
acoustic waves undisturbed by boundary reflections. This implies tanks 
of dimensions greater than ten wavelengths (about 15 meters at 1 kilo- 
hertz) for underwater testing. Coupler reciprocity is limited by the 
requirement to encase the entire active surface of two test transducers 
in a rigid-walled cell of dimensions less than 1/10 of a wavelength; 
useful underwater low frequency transducers limited to such dimensions 
can be difficult to build. 

The practice of reciprocity cal.ibration is not, however, 
constrained to the two conventional geometries. Rudnick [Ref. 7] 
discusses three compact geometries with known reciprocity factors; 
Bobber's reference text on measurements [Ref. 8] describes several 
more. These unusual test environments extend the practice of reci- 
procity calibration beyond the limits of of the two conventional geome- 
tries. A summary of these methods is included in Table I; each of the 
reciprocity parameters listed is the quotient of an area characteristic of 
the system geometry, and of the specific acoustic impedance of the 
medium . 

The compliant tube test system described in this paper adds 
another option to the array already available to the experimenter. It 
has particular value to the underwater acoustician because it 
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TABLE 1 

Some Reciprocity Calibration Systems 



Geometry 



Reciprocity Parameter 



Reference 



Free field 

Close— coupled 

Plane— wave 

Cylindrical— wave 

Helmholtz resonator 

Waveguide, 

resonant plane wave 



2dX 

PqC 

2TrALf 

PqC 

2A 

PqC 

2LV5a 

PqC 

27rALf 

Poc'Q 

ttA Lf 

PqC^ 



4 

5 
8 
8 
7 
7 



significantly lowers the frequency limit for underwater reciprocity cali- 
bration below the practical limits of the conventional geometries. 
Unlike rigid-tube methods, it can be accomplished within a convenient, 
inexpensive test cell, and although still limited by the requirement for 
accurate measurement of the system quality factor, does escape resonant 
reciprocity calibration's usual restriction to discrete frequencies. The 
'Virtual rigid walls" created by the standing wave pressure antinodes, 
and in which the transducers are placed for the test, can be adjusted 
to any resonant length and frequency by simply raising or lowering the 
water level. 

Our development of the compliant tube reciprocity parameter 
must begin from an understanding of wave motion in the water column. 
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and of the oscillatory expansion of the tube's circumference. From 
these equations, we will develop the energy balance of the oscillating 
system; this in turn will provide the result for the reciprocity 
parameter. 

8. COMPLIANT TUBE WAVE EQUATION 
1 . Waves in the Water 

Lighthill's development [Ref. 9] of the wave equation for a 
lossless fluid moving coaxially in a lossless tube of varying cross 
section includes the case of the compliant elastic tube in which our 
experiment is performed. The momentum and mass conservation equa- 
tions for an inviscid, ideal fluid are 



\dt dz/ dz 



( 1 . 6 ) 



d(po) 

d t 




= 0 



(1.7) 



where p is the acoustic excess pressure, u the acoustic coaxial particle 
velocity, and a the local time-varying cross-sectional area of the tube, 
which varies slightly about its equilibrium value A. When fluid entropy 
is constant, both the local density p and the local cross section a may 
be expressed as functions of the acoustic pressure. (The latter relation 
is just the inverse of the function which relates pressure to cross 
section variations.) Discarding terms of second order and higher 
reduces these equations to 



(1.8) 

^ at az 

a(po) . au (1-9) 

— — = -poA — 
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Taking the spatial derivative of (1.8) and the time derivative of (1.9), 
combining the results with the linearized equation of state p = (p -p(j)c^ 
and discarding terms above first order gives a wave equation 



_ -2 9^p 

dz'~ ^ 



where the wave speed c is defined by 



( 1 . 10 ) 



= Po 



1 d(po) \ _ 



PoA d p 



Po 



L '^p 

Podp. 



L — 

A 3pJ 



= Po K + D 



( 1 . 11 ) 



(The subscript "0" above implies evaluation of the derivative at the 
undeformed wall position and under adiabatic conditions.) The first 
term of equation (1.11) is the adiabatic compressibility or reciprocal 

bulk modulus K; if there were no dependence on cross section, this 

would reduce to the equation for the sound speed in the bulk medium, 
if the second term, called the distensibility D, is much larger than the 
adiabatic compressibility, the speed of sound waves in the tube will be 
determined by the characteristics of the wall rather than by those of 
the bulk fluid, and will be much slower. This in turn will reduce the 
size of the resonant tube required to operate at a given frequency. 

The harmonic solution of equation (1.10) with a pressure 

release boundary at z=0 and a rigid boundary at z=L, and for frequen- 

cies below the cutoffs for azimuthal and radial modes, is 



p(z,t) = P sin(kz)sin(cjt) 



( 1 . 12 ) 



with a corresponding velocity field 
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0.13) 



u(z,t) = cos(kz) cos 

PoC 



(cjt) 



The wave number k is limited by the boundary conditions to the 
discrete values 



k 



2tt _ 7r(n + 1/2) 
X ~ L 



n 0 f 1 f 2 f 3 f • • • • 



(1.14) 



This sound field has associated with it both kinetic and potential energy 
density. The kinetic energy of axial water motion in a thin radial cross 
section of length dz is 



dE«, 



^Adz 

2 



(1.15) 



The potential energy in the same section is found by considering the 
work done by the fluid as the tube boundaries expand: 



dEp 



r<r 


K ( 


\ .Ida 


p dV = 




Adz 


V A 


J n ^ 


' A dp 



A dz 

Po^pJ 



dp 



(1.16) 



= (K - Adz 



(1.17) 



where C is a coordinate of the periodic motion. 

2. Motion of the Wall 

The distensibility can be calculated from the material proper- 
ties of the tube wall and its geometry. In the case of a cylindrical tube 
which responds elastically to small excess stresses, the circumferential 
tensioning force component due to the pressure on a thin radial cross 
section of length dz is 
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(1.18) 



p/r/2 

F = 2Rdz p cos 0 d0 = 4Rdz p 

-n/3 

^ ^ A ~ 2hYdz ” hY ^ 



(1.19) 



This tensioning force causes a strain s defined by equation (1.19). 
Where p is small and a varies only slightly from. its undisturbed value 
A, the distensibility is then 




2R 

hY 



( 1 . 20 ) 



This stretching of the wail along its circumference has the characteris- 
tics of a damped, stiffness-controlled driven oscillator (neglecting long- 
itudinal effects). The ratio of the driving pressure at any cross 
section and the circumferential strain velocity is the mechanical impe- 
dance of the cross-sectional ring treated as a damped oscillator of 
negligible mass with elastic constant Y and damping coefficient R(m) 
[Ref. 10;pp. 10-12] 



2P sin(kz)sin((i;t— 5) 

£ = ^ n 9V) 

W = ^Rm+(Y/y)^ 0 = arctan^— (1-22) 

The polyvinyl chloride (PVC) tube used in this experiment is actually 
not elastic over any finite range of distortion, so the parameter Y is 
not Young's modulus; it is only a linear approximation to the plastic 
deformation that the PVC undergoes in harmonic motion. Within this 
linearization, the potential energy stored in a thin radial section of the 
wall as the tube is expanded by a small amount from its initial cross 
section can be found from the force-displacement integral 
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0 



F d£ = Dp^ Adz 

2 



0.23) 



with the force given by equation (1.18). 

C. COMPLIANT TUBE RECIPROCITY PARAMETER 

The arguments used by Garrett, et £[. , [Ref. 11] to derive the 
reciprocity parameter J for a rigid tube will also yield an expression for 
the reciprocity parameter of the free surface, compliant wall tube. As 
with the rigid tube, the solution depends on relating the energy 
content of the system E and the energy loss per cycle AE to the 
acoustic transfer admittance (reciprocity parameter) described by equa- 
tion (1.2). When the ratio of energy content to energy loss rate is 
known in terms of J and of primary measurements, the resonance 
quality factor equation (1.24) (which can also be found from primary 
measurements) can be used to evaluate J. 



1 . System Energy Content and Losses 

The compliant tube system stores energy both in the water 
column and in the tube wall. Kinetic energy of motion of the water 
pressurizes the water locally, which in turn stresses and expands the 
tube diameter. The energy balance exchanges between the three 
degrees of freedom are governed by the work equations for pressuriza- 
tion (equation (1.17)) and for expansion (equation (1.23)), and by the 
kinetic energy equation (1.15). Summing these equations gives the 
total energy in a thin cross section 



f E 



(1.24) 



fh - ^ AE/277 
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dE = dEu + dEp + dE« = A dz j (K - D)|-+ Dp 



0.25) 

0.26) 



where equation (1.11) has been used to replace K and D with p^c^. 
The sound field in the absence of attenuation is described by equations 
(1.12) and (1.13); substitution gives a form which can be integrated 
over the length of the water column 





cos^(kz)cos^(cjt’) 



p2 ^ ^ 

H 2 sin‘(kz) sin‘(cjf)> dz 

2poC ) 



0.27) 



ALP- (1.28) 

~ 2 

4-poC 



This is identical to the result for a rigid resonant tube [Ref. 11]. 

In the steady state, the cyclic energy loss equals the work 
done by the driving transducer. If this transducer is located at a 
pressure antinode, pressure and velocity must be in phase at the 
transducer face ( since the impedance to flow at a velocity node is infi- 
nite). The energy, transferred to the system is the work done, which 
can be written as 



dZ 

at 




u® n ds 



= P U sin^(cjt) 



0.29) 



where U is the volume velocity amplitude of the transducer surface. 
Integrating over one cycle gives 




0.30) 
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2. Reciprocity Parameter 



Our calculations above have yielded both the system energy 
content and the system energy loss (actually, its throughput at steady 
state). The energy equation for Q (1.24) immediately produces 

Q_ __E _ P ALf 
2 tt AE U2poc"^ 

The ratio of volume velocity amplitude to pressure amplitude that 
appfears in this result is the acoustic transfer admittance between pres- 
sure maxima (equation (1.2)), and must also be the reciprocity param- 
eter for hydrophones positioned at separate maxima when the column is 
driven at resonance. Equation (1.32) follows; again, it is the same as 
the reciprocity parameter for the rigid tube. 



_ -n-ALf 

tl.32) 

The solution of equation (1.32) requires the experimental measurement 
of three resonance parameters: the resonance frequency f; the reso- 

nance quality factor Q; and the wavelength' of the resonance standing 
wave, in order to compute the sound speed C. The accuracy of the 
value computed for J depends on our ability to measure these parame- 
ters accurately. The methods for doing this will be discussed in the 
next chapter. 

In using equation (1.5) with the result derived above, we 
must ensure that the pressures and voltages implied are indeed the 
maximum amplitude values at the antinodes. The test transducers' 
output will be proportional to the average sound field over their active 
surfaces, which at its maximum will be less than that which would be 
generated by a uniform field at the maximum pressure amplitude P by 
the factor 
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(1.33) 



_ 1 X f ^ J sin(kL,/2) 

r,--5)cosk(z-z„)ds 

Measured voltage outputs must be corrected by this factor before the 
results are used in the calibration equation when transducer dimensions 
are an appreciable fraction of a wavelength. For the LC-lO’s used in 
our experiments, this factor is 0.05 dB at 1kHz; for the LC-32's, it is 
0.27 dB. 

3 . Limits of Derivation 

The derivation of the compliant tube reciprocity parameter has 
been done under four important simplifying assumptions. To the degree 
that these approximations do not represent the real system, the results 
obtained by absolute calibration will be in error. 

The first of these approximations is the use of the solutions 
to the lossless wave equations in computing the total energy in equation 
(1.28). The PVC wall does attenuate the sound wave, resulting in a 
resonance sound field which resembles that derived by Kinsler, Frey, 
Coppens, and Sanders [Ref. 10: pp . 206-209] . Figure 1.1 shows the 

measured variation of the pressure amplitude with depth for the n=4 
mode. The pressure amplitudes at separate maxima will therefore not 
be identical, as is required by our development, which assumes an infi- 
nite standing wave ratio. 

A second approximation is the assumption of linear behavior 
for the polyvinyl chloride wall. In fact, PVC has no elastic region for 
strain, as was discussed in the derivation of the wall motion equation. 
The nonlinearity of the wall permits mode coupling into other resonance 
modes, particularly when two resonances have spatial pressure maxima 
near the same position and are related subharmonically . The result of 
such coupling when the column is driven at one frequency is to depress 
the amplitude of individual pressure maxima by an amount proportional 
to the degree of coupling at that position; again violating the require- 
ment for equal pressure amplitudes at the two test positions. 
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Third, the reciprocity concept requires that the receiving 
transducer approximate an infinite mechanical impedance to the medium. 
While piezoceramics do meet this approximation very well, the soft 
mounting material used the the hydrophone body and gas bubbles on 
the hydrophone surface do not. Their presence could significantly 
distort the sound field in the vicinity of the transducer. 

Finally, the projector used to drive the column into resonance 
does not produce a sound field identical to that generated by a small 
piezoelectric transducer suspended in the water column. The local 
sound field near the cylinder is some combination of radial modes which 
match the cylindrical boundary conditions to the overall plane wave, 
coupling the driving transducer to the standing wave system. The 
effects of this coupling factor are slightly frequency-dependent, and 
are not accounted for in this model. 

Our measurements demonstrate that the simple model derived, 
despite these deficiencies, produces consistent results to within ±0.8 
dB, but does not explain differences from a standard comparison cali- 
bration of up to 2.1 dB. The root-mean-square average difference is 
1.5±0.4 dB for the range of hydrophone types and frequencies tested. 
A method and apparatus for performing this type of calibration experi- 
ment, and the results of our experiments, follow. 
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Figure 1.1 Pressure Amplitude Depth Profile (n=4) 



II. COMPLIANT TUBE RECIPROCITY EXPERIMENTS 



A. RESONANCE PARAMETER MEASUREMENT 

The formula for the reciprocity parameter derived in the preceding 
section requires the accurate measurement of three resonance parame- 
ters: the resonance frequency, the sound speed at that frequency, 

and the resonance quality factor; Each of these measurements involves 
finding an extremum from a set of .data samples, which will inevitably 
contain noise. The experimenter who attempts to resolve extrema from a 
discrete sample set by comparing measured signal outputs performs the 
discrete equivalent of differentiation, which enhances noise and 
increases net error in the result. 

The solution to accurate, reproducible resonance parameter meas- 
urement which D. Conte and S. Garrett originally developed [Ref. 12] 
combines the high data rate of automated measurements with the power 
of curve fitting. The intuitive explanation of the effectiveness of this 
method is the human analogue; rather than take discrete measurements, 
experimenters instead "swing" the system slowly through the 

extremum, and let the eye interpolate the location of that point. The 
optimum automated processing analogue for this method is curve fitting 
to the minimum expansion of any function y(x) which will resolve 
extrema, the second order equation 



y(x) = Co -»• C,x Cjx^ 



The three coefficients of this curve can be found for any set of three 
or more measurements; if more than three measurements are made, the 
nonlinear regression methods described by Bevington [Ref. 13] can be 
used to reduce the noise error in the measurement of the extremum, as 
long as the span of the sample set around the extremum is small enough 
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that equation (2.1) is a valid approximation. The regression equations 
for the three coefficients are algebraic, and are derived in Appendix B 
(equations (B.12), (B.13), and (B.14)). The location of the extremum 
(again from Appendix B) is displaced from the coordinate origin of the 
sample set by 



Xm - Xo = -C1/2C2 



( 2 . 2 ) 



In pur application of this method to measurement of the resonance peak 
frequency , the sample set is composed of seven hydrophone output 
voltage amplitudes measured at equal frequency intervals; the sample 
set origin is shifted until the estimated peak frequency falls within one 
interval of the center frequency of the sample set. This is the peak 
frequency locating algorithm. Figure 2.1 is a example of this technique; 
the asterisks mark data points used to establish the coefficients for the 
curve and to compute the location of the maximum, and points marked 
with circles are discarded from the final fit because they are more than 
three steps from the extremum. The sound speed is found by 
sampling the pressure amplitude at regular intervals in depth from the 
surface, and again shifting the set used for curve fit until the esti- 
mated maximum or minimum is within one step of the center of the set. 
The positions thus located occur at quarter wavelength steps, and can 
be fit to 



Zm = rn\/4 + zo (2.3) 

to obtain the wavelength; the sound speed is then just the product of 
the resonance peak frequency and the wavelength. This is the sound 
speed measurement algorithm. 

The curve fit in frequency also yields a reduced-noise value for 
the peak amplitude at resonance, which is essential to accurate 



24 



measurement of the half-power frequencies required for the Q measure- 
ment algorithm (equation (1.24)). Once the half-power amplitude is 
known, the half-power frequencies can be found by sampling outward 
from the peak frequency until a pair of samples brackets this value at 
the minimum frequency interval permitted by the equipment (in our 
case, 1 Hz). The half power frequency is then calculated by linear 
interpolation . 

The curve-fitting method can resolve extrema in noise to a degree 
far surpassing the eye integration method it was drawn from. It has 
the additional advantages of sound mathematical basis, speed, reprodu- 
cibility, and independence from observer biases. It also has the consid- 
erable advantage of being easier on the observer. The method's main 
drawback is its bias toward finding the parabolic curve in any data, 
including pure noise; care must be taken in selecting stepsizes to 
ensure that the changes in the signal are not so small that they are 
overwhelmed by noise, and to ensure that the sample set span is not 
larger than the region about the extremum where a second order 
approximation is valid. The best method for checking new cases is 

graphical display of the data set and the curve for the operator 
approval . 

B. DESCRIPTION OF EQUIPMENT 

1. Control System ( Figure 2.3) 

The experiment is controlled by a venerable Hewlett Packard 
HP-85 laboratory computer with 32K memory extension, IEEE-488 bus 
controller, and appropriate ROM extensions. The program used is 
documented in Appendix C. Commands generated by the program over 
the data bus control the experiment's configuration changes during 
measurements through the two dual remote coaxial switches (Hewlett 
Packard HP 59307) and a Hewlett Packard HP 3421 Data Acquisition 
Control Unit. The latter unit is configured with two actuator relays 
and eight measurement channels (option 020). The two actuator relays 
control the operation of the reversing DC manipulator drive motor. 
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Figure 2.1 Parabolic Curve Fit in Frequency 



26 




Figure 2.2 ACTRE Test System 
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Figure 2.3 Control System Schematic 
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Figure 2.4 Transmitter/Receiver Schematic 
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The manipulator's screw drive, which is described below, also 
drives a Hewlett Packard HED 680 position encoder that increments a 
Hewlett Packard HP 5016 universal counter in totalizing mode. The 
changes in this totalized count provide (via the data bus) the program 
with a measurement of actual manipulator motion. 

2. T ransmitter ( Figure 2.4) 

A Hewlett Packard HP 3314 function generator provides a 
stable bus-controlled sine input signal at two voltage/current levels, 
and a synchronization signal at the selected frequency and constant 
voltage. The high voltage output is amplified and used to drive the 
piezoceramic transducer at rms voltages of 10 to 20 volts AC. The 
low-voltage, high-current demands of the USRD Jll projector [Ref. 14] 
are met by low-gain power amplification of the function generator's low 
voltage output. The transmit remote coaxial switch (714)^ selects the 
signal source from these two options and feeds the signal to the 
reversible transducer switch (706A), to the Jll projector, or through a 
test attenuator which returns a known signal to the receiver input 
switch (706B). 

A Hewlett Packard HP3478 digital multimeter, configured as a 
remotely controlled 5 1/2 digit AG voltmeter, is connected to the 
transmit signal line to measure the transmitting voltage. 

3. Receiver ( Figu re 2.4) 

The test transducers used in these experiments are; two 
Celesco LC-10 piezoceramic hydrophones (Figure 2.5), selected for their 
small size and good transmitting response; and a Celesco LC-32 piezo- 
ceramic hydrophone (Figure 2.6) of larger size. The voltage sensitivi- 
ties of all three units were measured by comparison with a reference 
standard USRD H56 hydrophone [Ref. 14]. The results are listed in 
Table II. The LC-lO's are mounted on 0.375 inch thinwall copper 
tubing; the LC-32 is mounted on 0.4375 inch tubing. Signal cables are 



^The remote switches are identified here by their IEEE-488 bus 
address codes. 
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TABLE II 

Comparison Calibration Resul-ts 



Frequency 

THz) 

740 

860 

1000 



(All values 

LC-10#2166 
& 2m cable 

-206.5 

-206.8 

-206.9 



0.8 decibels) 

LC-10#2338 
& 9m cable 

-206.8 

-207.3 

-208.8 



LC-32#85 
& 11m cable 

-212.3 

-213.4 

-214.3 



routed up the interior of the tubing to minimize effects on the sound 
field. The two test transducers are suspended from the manipulator 
platform. Their output is returned from the top of the manipulator 
assembly by two of three coaxial cables of the same length and meas- 
ured capacitance. The third of these cables connects the test attenuator 
into the receiver input switch (706B), which makes all receive signal 
path wiring nearly equivalent. The reversible transducer switch (706A) 
selects either its receive or transmit lines; the receiver input switch 
(706B) selects the lockin analyzer input from the reversible transducer 
receive line, the test hydrophone, or the test attenuator. 

The EG&G PAR model 5204 lockin analyzer filters out the 
components of the received signal which are in-phase and in-quadrature 
with the reference (synchronization) signal.' Their Pythagorean sum is 
used to obtain the magnitude. The analyzer is typically operated with 
an integration time of 300 milliseconds, which corresponds to an equiva- 
lent noise bandwidth of 0.8 Hz; this causes the lockin analyzer to act 
as a very sharp filter synchronized to the transmitter signal. The 
analog output of the model 5204 is passed to a measurement channel of 
the HP 3421 for digitization. 

4. Mechanical System 

The tests are conducted in a 6-inch (154 mm) inside diameter 
Schedule 40 polyvinyl chloride tube mated to the active surface of a 
USRD Jll electrodynamic underwater projector [Ref. 14], The tube 
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DIMENSIONS IN INCHES 



Figure 2.5 Celesco LC-10 Hydrophone 
( ReprintecT with permission of Celesco Industries, Inc.) 

stands 1180 mm above the Jll's active surface, and is filled to a depth 
of about 1 meter; its wall, which is smooth on the interior except for 
small regular extrusion marks, is 8 mm thick. The configuration 
resembles the USRD G40 calibrator used for comparison calibrations 
[Ref. 14]. (This system was originally constructed by Mills and 
Garrett for their work [Ref. 15] with fiber optic gradient hydro- 
phones.) Equalization for the J11 driver is provided by air pressuriza- 
tion of the installed air bag, monitored with a water-filled U-tube 
manometer. 

The test transducers are suspended on thinwall tubing 
supported on a rigid table which is raised and lowered on two one-inch 
rod tracks by a motor-driven screw drive; maximum range of motion is 
800 mm. The drive controls were described in the paragraphs on the 
control system above. 
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Figure 2.6 Celesco LC'32 Hydrophone 
( Reprintecf with permission of Celesco Industries, Inc.) 

C. PROCEDURES FOR EXPERIMENT 
1 . Preparation 

Before the test transducers are, suspended in the water 
column, the experimenter must measure the reversible transducer’s 
input capacitance, so that the transducer current can be computed from 
the transmitter input voltage and drive frequency. The active surface 
lengths of both hydrophones should also be measured for use in 
computing the transducer amplitude averaging factor (equation (1.33)). 
The test transducers are then suspended from the manipulator, wetted 
with a detergent solution, immersed in the water column, and brushed 
to remove bubbles from the active surface. If necessary, the tube's 
inner wall should also be brushed free of bubbles. The water column 
depth is measured with a meter stick. 
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2. Mode Selection 



Two parameters are available to vary the resonant frequency 
of the water column at which the calibration is done: the water depth 

L and the mode number n (equation (2.4)). The procedure requires 
wavelengths much larger than the coaxial lengths of the transducers. 
In particular, the wavelength should be much longer than the coaxial 
length of the reversible transducer T to ensure that T's surface is 
entirely within the velocity node and feels little or no impressed pres- 
sure when T is transmitting. The wavelength should also be at least 
four times the length of the test hydrophone H. Practical experience 
has shown that calibration results are most accurate when T and H are 
positioned at two different pressure maxima of the n=3 or higher order 
resonance modes. 



L = ^(n + 1/2) = Y, 



(2,4) 



3- Resonance Parameter Measu rements 

After positioning the test transducers at the approximate loca- 
tions of the deepest and next deepest accessible spatial pressure maxima 
with similar pressure amplitudes, the operator uses the peak frequency 
measurement algorithm to find the resonance frequency f. The trans- 
ducers are then withdrawn above the first pressure minimum. 

One of the transducers is now used to measure the locations 
of maxima and minima in the standing wave pattern. Linear regression 
of this data gives the quarter wavelength, which in turn gives the 
sound speed. 

Finally, the transducers are precisely repositioned at the 
maxima where the resonance frequency measurement was made. A 
precise measurement of the peak frequency and of the Q factor are 
made, completing the set of measurements needed to compute the reci- 
procity parameter from equation (1.32). 
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4. Calibration Measurements 



Before making the final measurements, the experimenter 
switches the attenuated test signal into the receiver output to measure 
the net effects of cable losses, preamplifier gain drift, and deviations 
between the lockin analyzer’s and the digital voltmeter's measurements. 
This correction is applied along with the amplitude averaging correction 
factors from equation (1.33) to the calibration voltage measurements 
that follow. 

The transmitter is switched to the reversible transducer T, 
and the receiver to the test hydrophone H. The transducer-to- 
hydrophone output voltage E(TH) is measured, corrected, and 
recorded, as is the transmit voltage V(T). The latter is used with the 
capacitive admittance of T to compute the transmitter current l(T). 

The transmitter is then switched to the projector, and the 
drive voltage adjusted to give a sound field of the same amplitude at H 
as was used in the previous step. The projector-to-hydrophone output 
voltage is recorded (it is equal to E(TH), within the limits of the 
projector drive adjustment error). The final calibration measurement is 
made with the receiver switched to T with the projector drive 
unchanged; the result is the projector-to-transducer output voltage 
E(PT). 

These measurements complete the requirements to compute the 
test hydrophone voltage sensitivity from equation (1.5). The same set 
of measurements can be used to compute the reversible transducer's 
sensitivity using 



= J 



^PT^TH 

EpH^T 



(2.5) 



The results for M are in volts per pascal; the customary units quoted 
are decibels re IV/yPa, which is found by calculating 



M = 20 log M - 120 



( 2 . 6 ) 
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5. Self-consistency Check 



When both transducer can be used as reversible transducers, 
the results of the reciprocity calibration can be checked promptly by 
exchanging the transducers’ roles in the calibration measurements. 
Repeat measurement of the resonance parameters and transducer reposi- 
tioning is not required. 

Unless the two transducers have identical nonlinear responses, 
this second set of readings will confirm the accuracy of the assumption 
of reciprocity. It will not, however, measure the error in computing 
the^ reciprocity parameter J, which must be estimated independently. 

D. RESULTS OF EXPERIMENTS 
1 . Reciprocity T ests 

Tests were run on the three transducers in all pair combina- 
tions at the three accessible modes: n= 3, 4, and 5. Modes below n=3 
were inaccessible because the LC-10 transducers’ transmitting response 
was too small for reliable measurements; modes above n=5 were attenu- 
ated into the 'tails' of the lower resonances, which obscured the half- 
power points needed to measure Q. Comparison of a number of 
reciprocity calibration test results showed reproducible results to within 
0.8 decibels. However, the calibration results in Table IV disagree 
systematically with the values obtained for the same hydrophones by 
comparison with the USRD standard hydrophone (Table II) by 1.5 deci- 
bels (Table V). 

2, Error Analysis 

The self-consistency error quoted in Table IV was determined 
by reversing the electrical roles of the two transducers without 
changing their positions; the difference in computed results for the two 
measurements is the degree to which the system is not reciprocal. The 
estimated error is the combination of the self-consistency error with the 
error in measuring the reciprocity parameter, listed in Table III. 
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The primary source of self-consistency error is the difference 
in pressure amplitudes at the two maximum locations used, which we 
believe is caused by attenuation. By careful selection of modes and 
locations, we were able to reduce this error to 0.2 dB. Another lesser 
source of self-consistency error is the changing effect of the support 
tubing used for the transducers on the sound field as they are lowered 
into the water column. Insertion of the transducers and their supports 
into the water column displaces about 0.4 percent of the water column 
in volume, causing a similar change in the column depth and in the 
positions of maxima and minima. The change in the minima is the most 
significant error encountered in measuring the wavelength of the 

standing wave; because the cosine function is flat to first order, 
however, the effect on transducer positioning at maxima is negligibly 
small . 

Compliant materials associated with the signal cables and the 
transducers appear to degrade the reciprocal nature of the geometry. 
This error is reduced to very low values by selection of transducer 
location and by eliminating compliant materials (signal cables, for 

example) in contact with the water. 

Parameter measurement error is dominated by the apparent 
variation in the system quality factor over the length of the tube. An 
increase of 6-10 percent in the value of Q determined by the half-power 
frequency method was observed as the depth of the measurement was 
increased by one wavelength; our model for the compliant tube system 

does not account for such effects, and is therefore limited in the accu- 

racy with which it predicts the reciprocity parameter to about 0.8 dB. 
This is the largest recognized source of measurement error in the 
compliant tube reciprocity method. 

The systematically lower sensitivities measured by the reci- 
procity method relative to those measured by comparison with the H56 
standard are not yet understood. This discrepancy may rise from the 
inadequacies of our model for the system energy balance, from the 
method used to measure Q, or even from technical difficulties with the 
comparison calibration. 



37 



3. Measurement of Wall Motion 



Since sound speed can be measured directly, it is not neces- 
sary to attempt the difficult measurement of the actual elastic strain in 
the tube wall as part of the calibration process. However, the degree 
of correspondence of the actual strain behavior to the model derived in 
the first chapter is an important indication of the validity of our 
method . 

The circumferential wall strain was measured by attaching 
constantan .250 inch strain gages [Ref. 16] to the exterior surface of 
the, tube at the depths of some of the pressure maxima located by the 
test transducer. A second lockin analyzer used as the sensing element 
of a resistive DC bridge measured the magnitude and phase of the 
strain oscillations sensed by the strain gage; by comparing the output 
of the strain gage with the magnitude and phase of the pressure meas- 
ured by the transducer at that depth, we computed the complex impe- 
dance (W,0) of the wall ring (equation (1 .21 )) . The results in Table VI 
agree in general with the measured sound speed in the tube. Also, the 
indicated damping phase effects are reasonably small, justifying our 
approximation for the wall potential energy computation of equation 
(1.23). 
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TABLE III 

Resonance Parameter Measurements 
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RMS average total estimated error AM - t 

M = Hydrophone voltage sensitivity (dB ro IV/yPa) 
AMr = Self consistency error (dB) 
am = Total estimated error (dB) 



Variations in Calibration Results 
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725 3/4 (2. 7, -80) 2.6 330±30 365±6 
868 3/4 (2.6, -76) 2.5 330±30 355±3 
1065 5/4 (1.9, -82) 1.9 310±30 338±5 



I n_. UTILITY OF COMPLIANT TUBE CALIBRATION 

A. LIMITS OF COMPLIANT TUBE CALIBRATION 

The compliant tube reciprocity method can be extended to low 
frequencies by simply increasing the length of the water column 
enclosed in the PVC tube, or by decreasing the wall thickness to diam- 
eter ratio (equations (1.11) and (1.20)). The practical lower frequency 
limi.t is the length necessary to support the n=3 standing wave, which 
from equation (2.4) is 



f 

■ 4 L 



(3.1) 



The upper frequency limit depends both on the onset of azimuthal 
modes, and on the increasing attenuation of the sound wave by the PVC 
wall. The former constraint breaks the plane wave symmetry when the 
frequency is greater than (0.586»L/R) times the fundamental longitu- 
dinal mode [Ref. 17]. The attenuation problem causes excessive varia- 
tion of the measureable system quality factor with depth and reduces Q 
to the point where half power frequencies are obscured by adjacent 
resonances . 

Test transducer geometry is restricted to shapes that are small 
enough in cross section that reflections from their lower surface do not 
significantly affect the standing wave. The length of the reversible 
transducer, as discussed in the procedure description, should be as 
short as possible consistent with adequate transmitting response. (A 
flat disk transducer might be more effective than the cylindrical LC-10 
we used.) The test hydrophone can be much larger, if the transducer 
amplitude averaging correction is properly computed. 

Accuracy of the method is primarily limited by the effects of the 
wall attenuation. A derivation which accounted for these effects, and 
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for the lesser effects of disturbances to the standing wave caused by 
the test transducers and their support tubing, could conceivably 
achieve much lower calibration uncertainities , 

A significant advantage of the compliant tube method is the use of 
pressure antinodes as virtual rigid boundaries at which to locate the 
test transducers. Many of the difficulties associated with mounting a 
transducer in an real boundary which must have some compliance are 
•thus avoided. 

B. POTENTIAL APPLICATIONS 

Two applications for the compliant tube system on a practical scale 
have been suggested to the author: a field/small lab calibrator, and a 

large, very low frequency calibrator. 

The portability and flexibility of the test system used for 
compliant reciprocity recommends the method for general calibration use 
in small laboratories and in the field. The 1-meter PVC tube provides 
calibration capability in the 1 kHz range; a 2-meter, 10-inch diameter 
tube would be more versatile and could reach frequencies of 500 Hz. 
The manipulator assembly need not be as complicated or massive as the 
one used in our experiment; any system which can position the trans- 
ducers within a few millimeters is adequate. Because the method 
involves resonance measurement methods as’ well as the techniques of 
reciprocity calibration, it may also be of interest for advanced instruc- 
tional laboratories . 

On a larger scale, a 10-meter, 30-inch diameter tube, which would 
fit upright in a two-story open bay, could reach frequencies below 100 
Hz and calibrate units similar in size to the Navy's TR-155 transducer. 
If current interest in very low frequency underwater acoustics indicates 
a need for an alternative to the traveling-wave tube system described 
in [Ref. 8], compliant tube reciprocity may provide the answer at a 
very low cost. 

There is also no reason we are aware of that the benefits of 
compliant waveguide usage should be limited to resonant calibration. 
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Construction of a compliant traveling wave tube calibrator with a non- 
reflecting termination or of a Helmholtz resonant calibrator using two 
compliant tubes should also be practical. The Helmholtz calibrator 
concept in particular has the potential of extending the range of accu- 
rate absolute calibration for the underwater acoustician to extremely low 
frequencies . 
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APPENDIX A 
GLOSSARY 



Subscripted variables in the text are represented using paren 



theses. 


due to word processor limitations. 


A 


Cross sectional area of the water column at rest. 


A(T) 


Active surface area of the transducer. 


a 


Instantaneous cross sectional area of the water column. 


CVO 


Parabolic curve fit coefficients. 


c 


Sound speed in a compliant tube. 


D 


Distensibility of the boundary wall. 


E(TH) 


Transducer-driven , test hydrophone output voltage. 


E(PH) 


Projector-driven, test hydrophone output voltage. 


E(PT) 


Projector-driven, reversible transducer output voltage. 


f 


Resonance frequency. 


H 


The test hydrophone used only to receive during the red 
procity calibration. 


h 


Thickness of the tube wall. 


l(T) 


Transmitter current in reversible transducer T. 


J 


Reciprocity parameter, or acoustic transfer admittance 
between two points. 


K 


Adiabatic compressibility of water. 


k 


Wave number at resonance = 2tt/\. 


L 


Depth of the water column. 


L(T) 


Transducer's active length normal to the plane wave. 


M 


Hydrophone voltage sensitivity. 


m 


Number of quarter wavelengths below the water surface. 


n 


Resonance mode index. 


n 


Unit vector normal to surface element ds. 


P 


Acoustic wave excess pressure. 


P 


Pressure amplitude of the acoustic wave. 


Q 


Resonance quality factor. 


R 


Inner radius of the tube. 



44 



R(m) Damping coefficient for tube wall vibrations. 

S Transmitting current response, 

s Surface area element. 

T The test transducer used both as a projector and as a 

hydrophone during the reciprocity calibration. 

U Volume velocity of the acoustic wave, 

u Particle velocity of the acoustic wave. 

V Volume element. 

W Complex mechanical impedance magnitude of the tube wall. 

x(0) Center coordinate of the data set used for curve fitting. 

x(m) Coordinate of extremum of the curve. 

Y Elastic coefficient for tube wall expansion. 

y(i) Sample set to be curve fit. 

Z(E) Transducer blocked electrical impedance. 

Z(m) Transducer open-circuit mechanical impedance, 

z Depth from surface of water column. 

z(0) Depth measurement offset. 

z(m) Depth of the mth pressure maximum/minimum. 

5 Sampling interval for curve fit sample set y(i). 

£ Areal strain in the wall. See equation (1.19). 

2T Length averaging correction factor. 

9 Complex mechanical impedance phase angle of the tube wall. 

X Wavelength of the standing wave. 

P Instantaneous local water density, 

p Ambient water density. 

i Transformation factor. 

Sum of squared errors in a curve fit to data, 
w Angular frequency at resonance. 

^ Areal expansion coordinate. 
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APPENDIX B 

PARABOLIC SMOOTHING ALGORITHM 



The algebra of deriving formulas for the best choice of coefficients of a 
nonlinear equation, such that the sum of squared errors described by 



x^ = i:xf=E 

is minimized , is described by Bevington in [Ref. 13], The minimum 
for the function defined by equation (B.l) is defined by the selection 
of coefficients which solve 



y(xi) - yi 



(B.l) 




ac 



= 0 



(B.2) 



where C is any of the coefficients of the curve to be fit to the data. 

The simplest function that will predict maxima and minima is the 
quadratic equation 



y(x) = Co + CiX + CjX' 



(B.3) 



The single extremum of this curve occurs at 



^0 — C1/2C2 



(B.4) 



with respect to the coordinate origin located at x. Let the sample set 
be seven measurements of the function y sampled at regular intervals 5 
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yi - y(xo + i<5) ; i= -3, -2. -1, 0. +1, +2, +3. 



(B.5) 



Evaluating the best fit function’s partial derivatives with respect to the 
three parameters for such a sampled data set gives equations (B.6) 
thru (B.11); summing and cancelling yield equations (B.12), (B.13), 
and (B.14). These analytic results can then be used with equations 
(B.3) and (B.4) to obtain the location of the maximum and its ampli- 
tude, or the smoothed value for the amplitude at any location. 



X-i — 96^Cz~ 33Ct+ Co“ y-s 



(B.6) 



X~z — 4<5^C2~ 2(5Ci+ Co~ Y-z 



(B.7) 



X-i — 6^Cz~ <5Ci+ Co~ y~i 



(B.8) 



X 0 — Co“ y 0 



(B.9) 



X + i = (5^C2+ c5Ci+ Co- y+i 



X+2 = 46^C2+ 2<5C(+ Cq— y+2 
X + i “ C 2 "t 36C)+ Cq— y+i 



(B.IO) 



(B.n) 



843^C2 = 5(y.,j+ y.^) - 3(y*,+ y_,- 4yo ) 

28(5C, = 3(y+j- y.j) + 2(y.,2- y- 2 ) + (y+i- y-i) 

2 IC 0 = 7yo + 6(y^,+ y_,) + 3(y^2+y-2) - 2(y + j+ y_j) 



(B.12) 



(B.13) 



(B.14) 
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APPENDIX C 

ACTRE: HP-85 BASIC PROGRAM FOR RECIPROCITY EXPERIMENTS 



A. WHAT THE PROGRAM DOES 

The Automated, Compliant Tube Reciprocity Experiment (ACTRE) 
program provides the operator with a menu of soft key controlled 
"tools” for resonance parameter and calibration measurements. Its 
desjgn was based on the need for maximum flexibility, rather than 
maximum unattended operation. 

All of the measurement "tools" are subroutines called from the 
HP-85 function keys by the operator. The major procedures are: the 
peak frequency location algorithm (KEY 1); the spatial optimum location 
algorithm, which also computes wavelength and sound speed (KEY 2); 
the Q measurement algorithm (KEY 5); and the reciprocity calibration 
measurement sequence (KEY 7). Supporting procedures provide the 
actual control statements for remote switching, function generator 
setup, manipulator positioning of the hydrophones, receiver output 
measurement, and graphics output. Any procedure may be terminated 
by calling the process clear sequence (KEY 8). Comments on the 
detailed structure of the program are included with the source code. 

The program was developed to fit the demands of the experimental 
apparatus , the experimenter, and the research goals of this work as well’ 
as the limits of the host computer permitted. Time and core memory 
constraints severely limited efforts to make the program friendly or 
accessible to the uninformed. 

B. TEST EQUIPMENT BUS CONNECTIONS 

The code is designed to use the IEEE-488 bus-controlled test 
equipment described below. All of the primary measurements are 
made- with bus-controlled equipment, except for the hydrophone output 
voltage measurement by the EG&G 5204 lockin analyzer, whose output 
signal must be converted using the HP 3421 D/A unit. 
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Bus assignment numbers used in the program: 



705 HP7470 hardcopy plotter. 

706 HP59307A remote controlled switch: 

A switch: reversible transducer; 

B switch: receiver input. 

707 HP 3314 function generator. 

709 HP 3421 Data Acquisition/Control Unit, with Opt 020 

board in slot 0.- Channel assignments: 

0: Manipulator motor drive up actuator 
1: Manipulator motor drive down actuator 
5: 5204 lockin analyzer magnitude signal. 

714 HP 59307A remote controlled switch (transmit side). 

720 HP 5316 Frequency counter, in count totalize mode. 

723 HP 3478 digital multimeter, configured as an AC volt- 

meter. 



C. 



VARIABLE ASSIGNMENTS 
1 . Defined Functions 

FNH(XI) Length averaging correction function (equation 

(1.33)). 

FNK(KO) Cycles KO through the indices 1,2,3. Used in plot 
scaling algorithm. 

FNL(XO) Returns decibel value of XO = 20 log(XO). 

FNM(XO) Returns voltage sensitivity in dB re 1V/uPa (equa- 
tion (2.6)). 

FNR(XO) Returns value of XO rounded to three significant 
places; required for compatibility with HP3314 function 
generator amplitude command. 
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2. Variables 



A(7) Amplitude sampling array, used for curve fitting algo- 

rithms . 

A2 Amplitude ratio of hydrophone output voltage to 

transmitter input voltage. Measured by the Amplitude 
Read procedure. 

A4 Half power amplitude in Q measurement procedure. 

A5 Lockin analyzer output signal to D/A unit (709) 

channel 5; ranges from 0-1 V to indicate 0-full scale 
magnitude reading. 

A7 Resonance peak amplitude, computed from curve fit. 

A8 Preceding measurement storage--in Q measurement 

procedu re. 

B(3) Plot scaling parameter array 

B1-B9 Plot parameters 

C1-C5 Sums used in wavelength measurement least-squares 

procedure [Ref. 13]. 



Cl : 


Sum 


of 


N9. 


C2; 


Sum 


of 


N9*N9 


C3: 


Sum 


of 


L9. 


C4: 


Sum 


of 


L9*N9 


C5: 


Sum 


of 


L9*L9 , 



C6 The quantity A in the least-squares fit. 

C7 Resonance sound speed (m/s). 

C3 Error estimate for C7 (m/s). 

D1 Water column depth(mm) 

D2 Water column cross section (mm~) 

D3 Water column volume(m^) 

D5 Length of T (mm) . 

D6 Length of H (mm) . 
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El 

E2 

E3 



E(PT); Output voltage of T , P driven (V) . 

E(PH); Output voltage of H , P driven (V) . 

E(TH); Output voltage of H , T driven (V) . 



All frequencies are in hertz(Hz). 

FO Center frequency of sample set. 

FI First peak detection / Lower 1/2 power point. 

F2 Second peak detection / Upper 1/2 power point. 

F3 Peak frequency estimate. 

F4 Lower search boundary. 

F5 Upper search boundary. 

F6 Drive frequency; argument passed to Amplitude Read 

Setup procedure. 

F7 Resonance peak frequency. 

F8 Error estimate for F7. 

F9 Stepsize. 

G2 Receiver preamplifier gain. 

G5 Lockin analyzer output expand setting. 

G7 Calibration attenuator value.' 

G8 Drive boost factor for preamp calibration test in 

Calibration routine. 

HI Length averaging correction factor for T. 

H2 Length averaging correction factor for H. 
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10 General use 

13 Size of plot files X(), Y(). 

14 Selection index for Operator Input procedure. 

15 Number of time constants T5 to delay before measure- 
ments . 

18 Select index for Switching procedure. 

19 Wavelength measurements count. 

J3 Reversible transducer(T) transmit current (A). 

J7 Reciprocity factor (V*A/Pa^). 

KO Water density (kg/m^). 

K3 Plot scaling array index. 

L9 Pressure maximum/minimum position (m) . 

Ml Hydrophone voltage sensitivity of T (V/Pa). 

M2 Hydrophone voltage sensitivity of H (V/Pa). 

N0-N2 Position counts; used to compute manipulator motion. 

N9 Number of 1/4 wavelengths from surface corresponding 

to L9 . 

Q7 Measured quality factor of resonance. 

RO Position encoder count ratio for down motion 

(counts/ mm) . 

R1 Manipulator drive rate (counts/sec), down motion. 

R2 Manipulator drive rate (mm/sec), down motion. 

R3 Allowed projector drive adjust error. 
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R4 

R5 

R7 

R8 

The S 

50 

51 

52 

53 

54 

55 



56 

57 
S9 



Allowed peak frequency detection range. 
Allowed min/max detection range. 

Q measurement coarse stepsize factor. 
Manipulator up/down motion rate ratio. 

variables are all control flags/indices. 

General use. 

Read enable flag. 

Manipulator motion stop flag. 

Plot option enable flag. 

Fast display refresh flag. 

Drive voltage select index; 

0: use Jn projector drive voltage V4; 
1: use piezoceramic drive voltage V3. 
Automatic drive voltage control enable flag. 
Resonance peak detection index. 

Plot execute flag and abscissa select index. 



All times are in milliseconds. 

T2 Manipulator motion time adjust. 

T5 Lockin analyzer time constant. 

T7 Additional delay time before measurement. 

T8 Size of additional delay T7 to be used when required. 

UO Curve fit amplitude at center of sample set. 

U1 Curve fit slope*stepsize . 

U2 Curve fit curvature *(stepsize) ^ . 

U9 Stepsize used in last curve fit ( for plot procedure). 
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All voltages are in volts(V) rms, except where indicated 



VI 


Transmitter input voltage. 


V2 


Hydrophone output voltage. 


V3 


Piezoceramic drive voltage(Vp-p) . 


V4 


jn projector drive voltage(Vp-p) . 


V5 


Lockin analyzer fullscale setting; input as milli- 
volts, and converted to volts. 


V6 


Function generator drive voltage (Vp-p). 


V7 


Transmitter input voltage at resonance. 


V8 


Curve fit amplitude at F7, returned by 5-point 
Voltage Measurement procedure. 


VV0-W3 


General use working variables. 


X(50) 


Abscissa plot storage array (Hz or mm). 


Y(50) 


Ordinate plot storage array (mV rms) . 


XO 


Reference x-coordinate for curve fit plot (Hz or mm). 



All variables beginning with are in millimeters(mm) 



Z(7) 


Sampling positions for X Measurement procedure. 


ZO 


Manipulator table position. 


Z2 


■Estimated position of pressure minimum/maximum . 


Z4 


Position of center of sample data set/ operator-ordered 
position to move to. 


Z8 


Measured stepsize in X Measurement procedure. 


Z9 


Ordered stepsize. 



3. Character String Variables 



AS(128) 


Title; used to identify printout and plots. 


M$(20) 


Code identifier. 


SS(5) 


Remote switch position identifier. 


WS(32) 


General character input. 
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D. SOURCE CODE 



1000 I progra/n actre-7 

1010 ! Automated, Compliant Tube Reciprocity Experiment 
1020 ! M.B. Johnson, Dec 84 
1030 I 

1040 PRINTER IS 2 @ OPTION BASE 1 @ DEG 
1050 DIM A$C128] ,H-fC5] ,W$C32] ,SIC5] 

1060 A$-"ACTRE-7 mod 8 " @ CLEAR @ DISP AT 

1070 ! Functions 

1090 DEF FNH(Xl) = X 1 / S IN ( RTD ( X 1 ) ) 

1090 DEF FNK(Kl) = K1 MOD 3+1 
1100 DEF FNL(Xl) = 20+LGT(Xl) 

1110 DEF FNM(Xl) = FNL(Xl)-120 

1120 DEF FNR(Xl) I rounds to 3 significant figures for HP'ITH 
1130 Wl=.01 



1140 
1 150 
1160 
1170 
1 ISO 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 



! begin loop 

IF IP (XI ) < l' THEN X1 = X 1 + 1000 (? Wl-Wl/1000 0 GOTO 1140 
IF IP(X1)>9 THEN X1=X1/10 0 W1=W1+10 0 GOTO 1140 
! end loop 

FNR=Wl+IP(Xl+100+.5) 

FN END 



I end : f uncti ons 



! Variable initialization 



REAL A(7) ,A2,A4,A5,A7,A8,C1 ,C2,C3,C4,C5,C6,C7,C 
REAL D3,E1 ,£2,E3,F0,F1 ,F2,F3,F7,G2,J3, J4, J7.K0, 
REAL N0,N1 ,N2,U0,U1 ,U2,V1 ,V2,V8, Y3,Z (7) , 20,Z2, Z 




SHORT B(3) ,31 ,B2,33,34,B5,B6,F4,F5 
SHORT R0,R1 ,R2,R3,R4,R5,R7,R8,T2,T 



,F6,F8,F9,.G5,3-’,32 

5,T7,T3,U9,V3,V4,V3,y6. 



1270 SHORT X(50) ,X0,Y(50) ,29 



12S0 


INTEGER BS 


,B9,I0,I 


4 


H, 


15,13 


T 0 

1 ‘ ? 


k:,k 


n 
; , 


,N3,N4, 


N9 


1290 


INT 


EGE 


R SO 


,S1,S2,S 




S4 , 


S5 , So 


) J 


99 








1300 


' i 


r i t 


ial 


values 


















1 3 1 0 


31 , 


S2, 


36 = 1 




















1320 


C 


S4, 


S9=0 




















1 3 3 0 


D2^ 


PI + 


i54+154/4 ' 


















1340 


G2-'= 


1 0 


l: / -• 


2S9. 639 


0 


GS = 


30 












1350 


R0 = 


512 


. 5 @ 


R 1=600 


0 


R2 = 


. 6/RO 


a ^ 


5 ~ 1 


0 


T2=tS3 




1360 


R 3 “ 


.01 


' F 


adjust 


or 


r or 














13";' 


R4, 


R5 = 


1.01 


' p e a k 


de 


tec 


t ‘-2 f r 


or 










1380 


14 = 


6 0 


GC3 


UE 2060 


















1390 


I 4 


3 a 


zo= 


0 a GOSUB 


206 0 ‘ p c 3 > i 


: on 


i [ 


iput 




1 4 0 '1' 


14 = 


2 0 


G03UB 2060 


1 


dep 


t h , te 


in p e r 


at Ur- 


e 






1 4 1 0 


R7 = 


120 


' 0 


outstep 


rati 


0 












14 20 


S ( 1 


1 0 
' ) 


(21 = 


2 0 B(3) 


•Z 


cr 

» J 


! p 1 0 


t sc 


al i n 


9 


p a r- a m e 


t er s 


1430 


15 = 


5 0 


T8 = 


0 ! read 


d 


e 1 a y f a c 


tors 










1440 


K0 = 


998 


. 2 ' 


k g / m ■ 3 


water 


density 
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1450 ! begin: equipment setup=========-=~-========-==========" 

1460 OUTPUT 709 ;"RS" ! deer A/D unit 

1470 OUTPUT 707 ;"PR" ! clear 3314 Ten generator 

1480 OUTPUT 723 ; "H2" ! xmi t vol tmeter 

1490 OUTPUT 707 ; "APIOMV" 

1500 OUTPUT 709 ; " F 1 R0RA02 1 N5LS5 ; T 1" ! a/d receiver output 
1510 OUTPUT 720 |"FN12" ! position counter 
1520 OUTPUT 709 ; "CLS5" 

1530 14^1 @ G0SU9 2060 ! LIA parameters 

1540 ! end: equipment satup = ” = = = = = = = = = = = = = = = " = "'- = = -“" = " = "= = -==" 

1550 ! begin: Fresh slate for processing = = = = = =■. = = = = = = = = r. = -- r-rr 

1560 I3,I9,N9,N4 = 0 ! plot file J<;vl/4 file sices ,wl/4 count 

1570 C1,C2,C3,C4,C5,C6=0 ! wavelength sums 

1580 F7,F3,C7,C3,Q7, V6=0 

1590 Ml ,M2=1000000 1 0 dB re IV/uPa 

1600 V4=.0001 3 03=10 3 T7=0 

1610 86=1 ' ADC on 

1620 !S=1 3 GOSUB 8120 ! switch 

1630 GOSUB 1300 1 set main menu 

1640 ! end: fresh slate for proces3ing = = = = = = = = = = = = = = = = = = = = = = --'' 
1 650 ! begi n : process! ng stage = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = '--~ 
1660 ! hold here for operator inenu selection 
1670 CLEAR 

1630 DI3P USING 3590 ; AT , F7 , F3 , C7 , C8 , 07 , FNM (Ml ) , FNM (M2 ' 

1690 IF S3 THEN DISP "plot on ";EL3E DI3P "plot of^ 

1700 IF 36 THEN DISP "ADC on" ELSE DISP "ADC off" 

1710 DISP ST ; " 20=" | ZO 

1720 KEY LABEL 

1730 FOR 17=1 TO 49999 

1740 IF S4 THEN 1760 i fast refresh on return from key routine 
1750 NEXT 17 
1760 S4=0 
1770 GOTO 1660 

1730 ! end: processing c t a g e = = = •=•-= •• " = = •= = = = •• - = - “ - ■: •• = - ;• •• - ■• • • - 



\ 
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1790 

1800 

1810 

1820 

1830 

1840 

1850 

1860 

1870 

1380 

1990 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1920 

1990 

2000 

2010 

2020 

2030 

2C40 



! = = = = = = =: = = = |*|ENUS = = = =-- = = = == = = = = - = = = --=-’ = = 
! begin; main msnu======== "========== 

ON KEYS l,"Fn" GOSUB 2650 
ON KEY# 2,CHRf(ll) GOSUB 4150 
ON KEY# 3, "move" GOSUB 7990 
ON KEY# 4, "opr" GOSUB 2060 
ON KEY# 5,"Q" GOSUB 3380 
ON KEY# 6,"a/pt" GOSUB 8120 
ON KEY# 7, "CAL" GOSUB 5190 
ON KEY# 8, "clear" GOTO 1930 
S4=l ! fast refresh 
RETURN 

! end: main ,Tienu" = = = = = = = = = = = = = = = " = = = = = 

! begin: clear slate menu 

! provides ax it/refresh processing 

OFF TIMER# 1 0 OFF TIMER# 2 

OUTPUT 709 ;"0PN0" @ OUTPUT 70? ;"0PN 

S1,S2=1 I flag reset 

13, S3, 39 ~0 ! plot dump 

ON KEY# 1 , "main" GOSUB 1600 

ON KEY# 3, "wipe" GOTO 1550 

ON KEY# 6, "quit" GOTO 8E40 

PRINT "OPERATOR CLEAR/RESET" 

S4"l I fast refresh 

GOTO 1650 I enter processing stage 

! end: clear menu""" = = =^“-~~~~-~'-~~“~ 
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2050 I ======PR0CEDURES====================================-= 

2060 ! begin: operator input procedure (HP-HUMAN bus ) ======-= 

2070 IF 14=0 THEN 2030 ELSE 2130 
2080 ! menu display 

2090 CLEAR @ DISP USING 3600 @ DISP USING 3650 
2100 INPUT WO 

2110 IF W0>9 THEN Z0=W0 @ W0=0 ELSE I4=W0 

2120 IF 14=0 THEN 2630 ! return to main menu 

2130 ON 14 GOTO 2140,2200,2260,2310,2360,2400,2530,2590,2440 

2140 ! 5204 LIA set t i ngs ( i 4= 1 ) 

2150 DISP "LIA Tullscals(mV),tim9 constant! ms), expand"; 

2160 INPUT V5,T5,G5S V5-V5/ ( 1000»65) 

2170 DISP "delay: H time consts, extra time(ms)";a INPUT 13, TS 
21G0 PRINT USING 8610 j 1 000 + V5 , G5 , T5 , 1 5 , TS 
21°0 GOTO 2630 

2200 ! water depth , temp ( i 4=2) 

2210 DISP "column depth(!nm),temparature(C)";a INPUT D1,W1 

2220 D3 = DUD2/ 10^7 

2230 PRINT USING 3620 ; 31, D3 

2240 PRINT USING 3640 ; W1 

2230 GOTO 2630 

2260 ' drive voltages(i4--3) 

2270 DI3? "drive lor P/T " ; V4 ; " / " ; V3 ; 3 INPUT Wl 

2280 IF LEN(WT)>0 THEN V4 = VAL(WI) Q V3-VAL ( WT I PCS ( WT , " / " ) ' 1 I ; 

2290 PRINT "P dri ve: V4; "0" 8 PRINT "T dr i ve ; " ; 03 ; " V " 

2300 GOTO 2630 

2310 ' comments (14=4) 

2320 DISP "comments" 

2330 INPUT W$ 

2340 IF LEN(W$)>0 THEN PRINT "#";WT 8 GQ'O 233) ELSE 2350 
2350 GOTO 2630 

2360 ' various constants(I4=5) 

2370 CLEAR 8 DISP "enter new values, then CONTINUE" G PAUSE 
2380 COPY 
:3'?0 GOTO 2630 
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2400 ! title (14=6) 

2410 AJ=AfCl,18] 

2420 DISP "date";@ INPUT Wl@ Af = AIS<WI 3 GOSUB 8070 
2430 DISP "run #" @ INPUT Af = AJ?iW^ @ GOSUB 8070 
2440 ! transducer chars only(i4=9) 

2450 AI=AIC1,643 

2460 H$ = "T is:" 0 DISP H$;@ INPUT W$0 A# = AI8<H$?a'lf 0 GOSUB 8070 
2470 DISP "T active 1 ength (mm) , capaci tance (pF) " | 0 INPUT D5,Y3 
2480 HI = "H is:" 0 DISP Hl;0 INPUT W$0 AI = AI?<HS .?<WI 
2490 DISP "H active 1 ength (mm) " j 0 INPUT D6 
2500 PRINT USING 3630 ; AI,D5,Y3,D6 

2510 Y3 = Y3/10'''12 § D5 = D5/ 1000 0 D6 = D6/ 1000 ! SI conversions 
2520 GOTO 2630 

2530 I plot enab 1 e/di sab 1 e ( i 4 = 7 ) 

2540 33=N0T S3 0 13=0 

2550 IF S3 THEN 2560 ELSE 2630 

2560 PLOTTER IS 705 

2570 LIMIT 32,232,30,130 0 LOCATE 22,122,27,97 0 PEN 0 

2580 GOTO 2630 

2590 ' po3ition(:4=8) 

2600 DISP "posi ti on=" ; ZO; 0 INPUT Wf 
2610 IF LEN(Wf)>0 THEN Z0=VAL(WT) 

2620 GOTO 2630 

2630 14=0 0 34=1 0 RETURN 

2640 ! end: o p e r a t c r input p r o c a d a r e = = = = = = = = = = - = 
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2650 ! begin! peak frequency measurement procedure====“===- 

2660 PRINT USINS S460 ; Al I title 

2670 CLEAR @ DISP "peak frequency measurement" 

2680 DISP " esti mate , sear ch r ange , st epsi : e" ; 3 INPUT F7,W0,F9 
2690 F4=IP(F7-W0) 

2700 F5^=CEIL(F7 + N0) 

2710 ! documentation 

2720 PRINT USING 3470 ; T I ME , F7 , WO , F 9 

2730 I warmup reading 

2740 IF S5 THEN V6=V3 ELSE V6=V4 

2750 F6=F7 3 T7=TS 0 G0SU8 6990 ! Kmtr on ,long delay 
2760 GQSUB 7310 ! warmup read 

2770 IF S3 THEN S9=l 0 GQSUB 5770 ! begin plot of a vs f 

2730 ! begin :initialice sampling array 

2790 F6-F4-3itF9 0 GQSUB 6990 

2300 FOR 12=1 TO 7 

2310 GQSUB 7310 ' read amplitude A2 

2320 IF 12(7 THEN F6=F6+F9 3 GQSUB 6990 ! xmtr on next ^raq 
2830 A(I2)=A2 
2340 NEXT 12 

2850 I end :fill sample array 
2360 ' b eg i n : r es onan C 2 search 

2870 S7 , F 1 , F2 , F7 = 0 ! peak detect flag and detection storage 
2380 FOR F0^=F4 TO F5 STEP F9 

2890 I fO is the freq at the center of the sample data set 
2900 F6 = F0 + 4 + F9 3 GQSUB 6990 ' :<mtr on next read freq 
2910 ! curvature at fO 

2920 U2=(5*(A(7)+A(l))-34(A(5)+A(3))-4*A(4))/34 
2930 DISP USING 3430 ; F0,U2 
2940 IF U2<0 THEN 2950 ELSE 3120 
2950 ' slope at fO 

2960 Ul = (3MA (7) - Ad ) ) ■'■2» (A (6) - A (21) MA (5; - A(3) : ) / 23 
2970 ! peak frequency 
2930 F3.= F0-F9-U1/ (2-U2) 

2990 DISP USING 3490 ; F3 

30)0 IF AES (F3---0X"9*R4 THEN 3010 ELSE 3 120 
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30i0 ! peak detected at center of sample array 
■3020 ! peak amplitude 

3030 U0=(7*A(4)+6*(A(3)+A(3) >+3*(A(2)+l^(6) )-2*(A(7)+A( 1) ) )/21 
3040 A7=U0-U1»U1/(4»U2) 

3050 S7=S7+1 3 V7=01 

3060 IF S5 THEN V3=V6 ELSE V4=V6 

3070 IF S7=l THEN F1=F3 ELSE F2=F3 

3080 PRINT USING 8500 ; F3, 1000*A74V1 , 1000»U2 , 1000*V1 
3090 U9 = F9 @ !(0 = F0 
3100 GOTO 3160 

3110 ! end:peak detected block 
3120 I begin:peak not detected block 
3130 S7=0 I unset peak detect flag 
3140 ! end: peak not detected 
3150 i end: tests at FO 
3160 DISP 

3170 IF S7<2 THEN 3130 ELSE 3240 ! end search on two detects 
3130 ! push down sampling array 
3190 FQ.R 10^1 TO 6 
3200 A(I0)-~A(I0+1) 

3210 NEXT 10 

3220 60SUB 7310 S A(7)=A2 
3230 NEXT FO 
3240 ! end search 

3250 IF S3 THEN G03UB 5330 ' show plot 

3260 IF S7^:2 THEN COPY e GOTO 3270 ELSE 3310 

3270 IF Fl^O THEN PRINT "search failed" 0 GOTO 3350 ELSE 323v 

3280 I one good reading 

3290 F7^=Fi a F3--^F9 

3300 GOTO 3330 

3310 I two good r e ..i d i n g s 

3320 F7:= (F1+F2) /2 6 F3 =ABS (F2-F 1 ) ! average of two good rn.dmys 

3330 ' final results 

3340 PRINT USING 3510 ; F7,F3 

3350 S 4 ~ 1 ' fast r e f r e s h 

3360 RETURN 

3370 I end: center freq mmt pr oc edur e = 
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3380 ! begin: Q measurement procsdurs="============== 

3390 CLEAR a DISP "Q measurement" 

3400 DISP "center freq ="iF7 
3410 Q7=0 

3420 PRINT USING 8570 ; T I ME , r 7 , Z 0 , 3$ 

3430 GOSUB 4950 ! measure peak voltage v3 at 47 

3440 A7-V8/V1 

3450 ! set search limits 

3430 F9=CEIL(F7/R7) 

3470 F4=MAX (0,F7-20*F9) 

3480 F5~-F7 + 20 + F9 

3490 DISP "peak : 1000 + A7»V1 

3500 A4= A7* . 707 1 0373 1 1 85 ! 1/2 power amplitude 

3510 DISP "half p ower : " ; 1 000*A4»V 1 

3520 l'begin:search for lower 1/2 point 

3530 F3^!P(F7) 

3540 A9=A7 

3550 I begin :loop to find f3<fl 
3560 F3=F3-F9 

3570 IF F3>F4 THEN GOSUB 3990 ELSE 4090 ' ;;mtr on, or 

3530 30SU3 7310 I measure A 2 

3590 DISP "q ( A2- A4 ) » 1 000>0 1 

3300 ' check for downhill motion 

3310 IF A2>A8 THEN 3320 ELSE 3350 

3320 ! uphill; low 1/2 pt obscured 

3330 PRINT USING 8520 ; F3,A2,A9 0 COPY 

3340 GOTO 4110 

3350 ! d 0 w n h 1 1 1 - r e s e t c h e c k v a 1 
3330 AS=A2 

3370 IF A2<=A4 THEN 3380 ELSE 3550 
3330 I end loop 

3390 ' begin loop: bracket 1/2 point 

3700 IF F3:F7 THEN F3=F3+I S G03UB 3990 ELSE 4090 

3710 A3=A2 ' save last sample 

3720 GOSUB 7310 ! new sample a2 

3730 DISP "q ; (A2 -A4) * 1000*01 

3740 IF A2>=A4 THEN 3750 ELSE 3390 

3750 ! end loop. bracks ted "interpolate 

3730 F1^F3-(A2-A4) / (A2'A8) 

3770 PRINT USING 3530 ; FI 
3730 ' end search for FI 
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3790 ! begin: search for hi 1/2 power pt F2 
3800 F6=IP(F7) 

3310 A8=A7 ! check value 
3820 ! begin: loop for f6>f2 
3830 F6=F6+F9 

3840 IF F6<F5 THEN GOSUB 6990 ELSE 4090 

3850 GOSUB 7310 ! measure A2 

3360 DISP "q ( A2-A4) *1 OOO^V 1 

3870 1 check downhill motion 

3380 IF A2>A8 THEN 3390 ELSE 3920 

3890 ! uphill-hi 1/2 obscured 

3900 PRINT USING 3520 ; F6,A2,A3 S COPY 

3910 GOTO 4110 

3920 ! downhi 1 1 -reset check val 
3930 A8=A2 

3940 IF A2<=A4 THEN 3950 ELSE 3320 
3950 ! end loop 

3960 1 begin loop, bracket 1/2 point 

3970 IF F6>F7 THEN F6=F6-1 § GOSUB 6990 ELSE 4090 

3930 A3=A2 ! save last sample 

3990 GOSUB 7310 I new sample a2 

4000 DISP "q! (A2-A4)MOOO*Ol 

4010 IF A2>=--A4 THEN 4020 ELSE 3960 

4020 ! end loop. brack a ted-interpolate 

4030 F2=F6+(A2-A4) / (A2-A3) ' 

4040 PRINT USING 3540 ; F2 

4050 G7 = F7/ (F2-F1 ) ! quality factor 

4060 S4=l ! fast refresh 

4070 PRINT USING 5550 ; 07 

4030 34=1 ! fast refresh 

4090 RETURN 

4100 ! b a g i n : e r r 0 r p r o c e s s i n g 
4110 PRINT USING 8560 ; F6 
4120 COPY 
4130 GOTO 4080 

4140 ' end:Q measurement pr ocedure = = = = --= = = = = 
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4150 ! begin: wavelength measurement pr ocedure======= 

4160 CLEAR 0 DIS? "min/max positioner" 

4170 IF S5 THEN V6=03 ELSE V6=V4 ! select T or P driv 
4130 F6 = F7 0 T7 = T3 0 G0SU3 6990 ! xrntr on, long delay 
4190 GOSUB 7310 ! warmup read 

4200 IF S3 THEN S9=2 0 GOSUB 5770 ! plot start, a vs 

4210 DISP "position="; ZO 

4220 DIS? "stepsizs";0 INPUT Z9 

4230 PRINT USING 3720 ; TIME,Z0,Z9 

4240 ! beg i n : i n i t i al i Z9 sample array 

4250 FOR 12=1 TO 7 

4260 GOSUB 7310 ! read amplitude A2 

4270 Z4=Z0 0 GOSUB 7610 ! move down :9 mm 

4230 Z{I2)-=Z4 0 A(I2)=A2 

4290 NEXT 12 

4300 ' end: initialize array 
4-310 S7 = 0 

4320 ' begin: search loop 
•^330 Z4 = Z(4) ! array center 
4340 I p a r a b 0 1 i c tit 

4 350 U2=’(5* (A(7) +A ( 1 ) ) -3 M A ! 5 ) +A ! 3 ) ) -4 » A ( 4 ) ) / 54 
4360 U1=(3*(A(7) - Ad) ). + 2^(A(6)-A(2) ) + (A (5) -A (3) ) ) ,'23 
4370 U0=(7^A(4) ^-6# (A(5) tA( 3) ) +3^dA(6) <-A(2) ) -2^ ;A;7) ^ A 
4330 U9 = Za 0 X0 = Z4 0 A7-U0-U 1 *U 1 / ( 44U2 ) 

439C Z2 = Z4-Z3-*Ul,^(2tU2) ! predicted optimum 

4400 DIS? USING 3660 ; Z4,Z2,U1,U2 

4410 IF S7 OR ABS(Z2-Z4)<R5tZ9 THEN 4420 ELSE 4460 

4420 ! hi t-' document 

4430 IF S7 THEN 4550 ! end search 

4440 37=1 ' set search end 

4450 Z3=Z2 ' save peak estimate 

4460 ! continue search 

4470 ' push down sample array 

4430 FOR 10=1 TO 6 

4 490 Z(I0)=Z(I0 + 1) 0 A(I0)=A(I0^-1) 

4500 NEXT 10 

4510 GOSUB 7310 ' read a2(z0) 

4520 Z(7)=Z0 0 A(7)=A2 

4530 IF 57 THEN PRINT USING 8730 ; Z 2 , 1 000 > A7 1 , U 1 , 'J 
7610 

‘^540 30T0 4320 

4550 ! end: search loop 
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4560 I print Hnal sstiiTiate 

4570 PRINT USING 8730 ; 22 , 1 000»A7*V 1 , U 1 , U2 

4530 Z2«(Z2+23)/2 ! average 

4590 PRINT "average="; Z2; "mm" 

4600 IF S3 THEN GO3U0 5930 ! show plot 

4610 I begin data review 

4620 IF 36 THEN 4630 ELSE 4640 

4630 IF S5 THEN V3^V6 EL.SE V4 = V6 ! store ADC results 
4640 DISP CHRf(ll)|"/4 c oun t = " ; N9 + 2 ; 0 INPUT Wf 
4650 IF LEN(Wf)>0 THEN 4660 ELSE N7=N9+2 0 GOTO 4680 
4660 W0=VAL(WI) 

4670 IF W0>0 THEN N9=W0 ELSE 4910 
4630 ! begin :accept measurement 
4690 19=^19+1 

4700 W1 ,L9=Z2/1000 ! meters 
4710 ! begin: suras 
4720 C1-C1+N7 ' ';< 

4730 C2 = C2 + N9^N9 ' 

4740 C3^C3+L9 I -y 

4750 C4 = C4 + N9+L? ' "-xy 

4760 03-^05 hL 9»L9 ' -y''2 

4770 C6=--I9 + C2-C1^C1 I del factor 

4730 I end: sums 

4 790 ' b e g 1 n : c 0 m p u t e sound s p d 

4300 IF 06)0 THEN 4810 ELSE 4370 

4 BIO ' two points cr more 

4320 (I?»C4-C1 *C3) /C6 I slope 

4330 NO" ( C2 *C3 -C MC4 ) /C6 ' intercept 

4340 C7-4M<I1+F7 

4850 IF I9>2 THEN 4860 ELSE 4870 

4 360 C3^4 + 3QR ( ! Cd-" I WO ♦ WO-t N i + N 1 » 2 2 -2 'W0-» 23 - 2^W 1 ^24*- 
19-2) ^H?/C6! *F7 
4 2 “ 0 ! document 

4330 PRINT USING 3690 ; I 9 , 27 , 23 , 4000 ^ W 1 
4390 GOTO 4930 

“900 i end ; a c c e p t m e a s u r e ra t 
4910 ! begin: discard raeasurerat 
4920 PRINT "mrat discarded" 

•1 930 S4-1 0 RETURN 

49 40 ' end: wavelength in e a s u r e m e n t p r o c e d u e ' = 
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4950 ! begin: 5-pt peak voltage measurement prccecure 
4960 IF 35 THEN V6=V3 ELSE V6=V4 ! select T or P drive 
4970 F6=F7 @ T7=TB @ GOSUB 6990 ! xmtr on ,long delay 
4980 GOSUB 7310 I warmup read 

4990 ! beg i n : f 1 ve-p oi n t mmt of peak amplitude v2 
5000 S0=S6 Q S6=0 ! ADC off 
5010 F6=F7-2 @ GOSUB 6990 ! set xmtr Fn-2 
5020 GOSUB 7310 ! read v2(Fn-2);lst terra 
5030 F6=F6+1 @ GOSUB 6990 ! xmtr on Fn-l 
5040 08 = -(3?^V2/35) 1 first term 
5050 GOSUB 7310 ! read at Fn-l 
5060 F6^F6+1 0 GOSUB 6990 ! xmtr on Fn 
5070 V8=V8+l2*V2/35 ! 2d term 
5080 GOSUB 7310 ! read at Fn 
5090 F6=F6+l 0 GOSUB 6990 ! xmtr on Fn+1 
5100 Va=V8+17»V2/3S ! center term 
5110 GOSUB 7310 ' read at FnM 
3120 F6=F6+1 0 GOSUB 6990 ! xmtr on Fn^2 
5:30 08=08+12^02/35 ! 4th term 
5140 GOSUB 7310 ! read at Fnr2 
5150 VS=V8-3»V2/35 ' result at f7 
5160 36= SO ! restore ADC 
5170 RETURN 

5130 ' end: peak voltage ni m t = = = = = •■= = = = = = = = = = =■• = = = = = -■ = ■=•- = =- = = -= 
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5190 ! begin: calibrate hydrophone procedurs = = = = -~ = = = = = = '=- = 
5200 W1=C7/F7 ! wavelength 

5210 H1=FNH(PI»D5/W1) ! length average corrn for T 
5220 H2 = FNH (P I if D6/W 1 ) I length average corrn for H 
5230 PRINT USING 8700 | A# , FNL (H 1 ) , FNL (H2) 

5240 I beginipreamp calibration 
5250 CLEAR @ DISP "calibration" 

5260 19=4 @ GOSUB 8120 ! switch 

5270 DISP "preamp calibrate" @ WAIT T3»3 

5230 F6 = F7 0 V6=V5#G3/32, 0 T7 = T9 0 GOSUB 6990 ; xnitr on 

5290 GOSUB 7310 ! read LIA out 

5300 G2=G7»A2 

5310 PRINT USING 3740 ; FNL ( G2 ) , 1 000 + V I 

5320 ! end: preamp calibration 

5330 ! T'H measurements 

5340 18=3 0 GOSUB 8120 ! switch 

5350 DISP "T}H" 0 INPUT Wl 

5360 GOSUB 4950 ' 5-pt mmt of v8 

5370 E3=y8*H2/G2 ! Eth 

5380 J3=Vl+2+PI+F7+Y3 I amperes 

5390 PRINT USING 3750 ; VI , 1000+J3 , 1000000 ^E3 

5400 ' end TIH 

5410 ' begin PIH 

5420 I8=-2 0 GOSUB 3120 

5430 DISP "PIH" 0 INPUT W$ 

5440 ' P drive adjust 

5450 F6 = F7 0 V6 = V4 0 T7 = T3 0 GOSUB 6990 : :<mtr on 
5460 GOSUB 7310 ! warmup reading 
5470 I begin: adjust loop 
5490 GOSUB 6990 ! xrntr on 

5490 GOSUB 7310 0 E2=V2+H2/G2 ' road v2, correct for preamp 
5500 DISP USING 3710 ; 1000+V6 , 1000 +V2/G2 , 1 000000 ' E3 
5510 IF ABS (E2-E3) :>R3 ‘£3 THEM V6 = r"NR ( V6 ♦E3/E2) 0 GOTO 54-;> 
5520 ' end adjust loop 
5530 V4=V6 

5540 PRINT "P drive; ";V4; "V" 

5550 ! end : F drive adjust 

3560 GOSUB 4950 ' vS mmt 

5570 E2=VG+H2/G2 ' Eph 

5380 PRINT " Eph = " ; r000000»£2: "uV" 

5590 ! end PIH 
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5600 ! begin P>T 

5610 DIS? "P}T" 0 I3=-l § G0SU8 8120 ! 

5620 GOSUB 4950 ! v8 mmt 

3630 E1=VS<H1/G2 ! Ept 

5640 PRINT " Ept = " ; 1 000000»E I ; " uV " 

5650 ! end P>T ramts 

5660 ! begin sensitivity calc 

5670 J7=PI»D3*F7/(K0^C7*C7»Q7) 

5630 PRINT USING 8760 ; F7,C7,Q7 
5690 PRINT USING 0770 ; D3,K0,J7 
5700 ^U=SQR(J7^^EUE3/(E2^^J3) ) 

5710 t12 = SQR(J7»E2*E3/(El*J3)) 

5720 PRINT USING 0780 ; FNM ( M 1 ) , FNM ( M2 
5730 ' end sensitivity comp 
5740 S4=l s RETURN 

5750 ' end:calibrate phones = = •- = = = = = = "• = 



switch 
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5760 ! ========PL0T PR0CS======= 

5770 ! begin: plot initialize procedure = = = = = = = -- = = 
5780 DISP "load plotter";® INPUT Wl ! dummy input 
5790 IF LEN(WI)>0 THEN S3=0 @ RETURN ! abort plot 
5800 PEN 1 S FRAME 0 PENUP 
5810 RETURN 

5820 ! end: plot initialize procedure============ 

5830 ! begin: show plot procedure" === === = ==’= = = = = 

5840 ! begin: scaling 

5850 81,84=99999 0 82,85=-Bl 

5860 FOR 10=1 TO 13 

5870 W0=X(I0) 0 W1=Y(I0) 

5830 IF W0>82 THEN B2=W0 
5890 IF W0<81 THEN B1=W0 
5900 IF HI >85 THEN B5=W1 
5910 IF WKB4 THEN B4 = N1 
5920 NEXT 10 

5930 DISP "data range x: ";B1;B2 
5940 DISP " y: " ; 84; 85 

5950 ! X axis scale 
5960 83=5 ! K-interval 
5970 WO =5*83 ! major divisions 
5980 Be=(CEIL(B2/W0)-IP(Bl/W0) ) +5 
5990 IF BB>50 THEN 6010 ELSE 6000 
6000 IF 88<20 THEN 6080 ELSE 6160 
6010 ! too many divisions 
6020 K3=l 

6030 B3=B3*3(K3) 0 K3=FNK(K3) 

6040 W0=5*B3 

6050 B8=(CErL(32/W0)-IP(Bl/W0) ) *5 
6060 IF B8>50 THEN 6030 ELSE 6160 
6070 I end too many divisions 
6080 ' begin too -few divisions 
6090 K3=: 

6100 83 = 83/8 (1(3) 0 K3 = FNi (K3) 

6110 IF 83)0 THEN 6120 ELSE 6200 
6120 W0=3*S3 

6130 SS=(CEIL (B2/W0) -IP (81/WO) ) *5 
6140 IF 89<20 THEN 6100 ELSE 6160 
6150 ! end too -few divisions 
6160 ' end if 

6 170 31 = IP(31/W0) »W0 I l ow :< bdry 
6180 82 = 81 + 83 '+38 ! hi x bdry 
6190 SOTO 6250 

6200 ! method failod-set default ;< scaling 
6210 B3=(B2-81) /lO 
6220 88=10 
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6230 DISP "default x -scaling" 

6240 ! end default scaling 
6250 ! end x scaling 
6260 ! begin y scaling 
6270 86^.002 ! nV 
6280 W2=5*B6 

6290 B9=(CEIL(B5/W2)-IP(B4/W2) )*5 
6300 IF B9>40 THEN 6320 ELSE 6310 
6310 IF B9<10 THEN 6380 ELSE 6450 
6320 ! begin too many divisions 
6330 K3=3 

6340 B6=B64B(K3) 0 W2=5*B6 § K3=FNK(K3) 
6350 B9=(CEIL(B5/W2)-IP(B4/W2) )*5 
6360 IF B9>40 THEN 6340 ELSE 6450 
6370 ! end too many divisions 
6380 ! begin too few divisions 
6390 K3=l 

6400 36=B6/3(K3) § N2=5+B6 0 K3^FNK(K3) 

6410 IF B6>0 THEN 6420 ELSE 6490 

6 420 B9=(CEIL(35/W2)-I?(B4/W2) )?(-5 

6430 IF B9<10 THEN 6400 ELSE 6450 

64'i0 1 end too few divisions 

6450 ! end if 

6460 B4--=-IP(B4/W2) M42 

6470 B5=B4+B6»39 I hi y bdry 

6400 GOTO 6540 

6490 ! m e t h 0 d ' f a i 1 “ 3 e t default y - s c a 1 i n g 
6300 D6-(B5-B4)/10 
6510 33=10 

6520 DISP "default y-scaling" 

6530 ' end default y scaling 

6540 ! end y-scaling 

6550 SCALE B1,B2,B4,BS 
6560 ' end scaling 

6570 ! plot points 

6530 OUTPUT 705 ;"SMo;" 

6590 K3=I3-7 

6600 IF K3:i THEN K3=0 0 GOTO 6640 

6610 FOR 10=1 TO K3 

6620 PLOT X(I0) ,Y(I0) ,0 

6630 NEXT 10 

6640 OUTPUT 705 ; "3M + ; " 

6650 K3 = K:3+1 

6660 FOR I0=K3 TO 13 

6670 PLOT X (10) ,Y(I0) ,0 

6680 NEXT 10 

6690 OUTPUT 70S ; "SM; " 

6700 CSIZE40FXDl,3 
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6710 LAXES B3,B6,B1 ,84,5,5 

6720 N4=N4+1 0 HI=AIC33,37] 0 PRINT "plot ";Hf;"#“;N4 
6730 SETGU 

6740 MOVE 0,62 0 LOIR 90 0 LORG 6 

6750 LABEL USING "K" ; " amp 1 i tude (mV rmsl" 

6760 LDIR 0 0 LORG 4 0 MOVE 72,12 

6770 IF S9 = 2 THEN Wl = "depth (mm) " ELSE = requency (Hz ) " 
6780 LABEL USING "K" ; W$ 

6790 MOVE 72,5 0 CSIZE 5 

6800 IF S9--=2 THEN LABEL USING 8670 ; F6,SJ,H$,N4 ELSE LAB 
8630 ; Z0,ST,H#,N4 
6310 SETUU 0 PENUP 
6820 1 end: point plot 
6830 ! begin: curve tit plot 
6840 B3=(B2“B1) 7100 
6350 FOR W1=B1 TO D2 STEP B3 
6860 W2=W1“X0 

6370 W3 = 1000*V1MU0 + U1*W2/U9 + U2» (W2»W2) / (U9-»U9) ) 

6880 IF W3>B4 AND W3<B5 THEN PLOT W1,W3 
6390 NEXT W1 
6900 FENUP 

6910 MOVE 0,0 0 PEN 0 
6920 ! end: curve tit plot 
6930 ! begin: clear plot 
6940 13=0 0 89=0 
6950lend;clearplot 
6960 RETURN 

6970 ! end: plot procedure============================== 
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6930 ! ======MEASUREMENT PR0CEDURES==================--======-- 

6990 ! begin: amplitude read setup procedure=====“===========“- 

7000 ! i nput : 36 , F6 , V6 out: timer SI set for read enable, read sn 
able si off 

7010 IF 36 THEN 7020 ELSE 7250 
7020 ! begin: ADC block 
7030 ! automatic drive adjust 

7040 OUTPUT 707 USING 8800 ; V6,F6 ! HP3314 drive on 
7050 S1=0 i ON TIMERS 1,T7+I5*T5 GOSUB 7550 i read disable 
7060 I begin: loop to adjust V6 until LIA output is in limits 
7070 WAIT 2*T5 

7080 ENTER 709 ; A5 ! read LIA magnitude-ignore read disable 
7090 IF A5>.8 THEN 7100 ELSE 7130 
7100 ! LIAout too high 

7110 IF ASM. 12 THEN V6 = FNR (V6* . 5) ELSE V6--FNR (V64 . 7) 

7120 GOTO 7170 

7130 IF A5<.1 THEN 7150 ELSE 7290 ! here's the loop exit point' 
7140 1 LIA out too low 
7150 V6=MIN(10,FNR(V6*15) ) 

7160 GOTO 7170 

7170 ! reset HP3314 and retest 

7180 OUTPUT 707 USING 8810 ; V6 

7190 OFF TIMERit 1 2 31=0 

7200 ON TIMERS 1 , T8+ ( 1 5 -2 ) *T5 GOSUB 7550 

7210 PRINT 3$n,l];" dr i ve" : V6; "V" 

7220 ! hold for read enable 
7230 IF SI THEN 7060 ELSE .7220 
7240 ! end: ADC block 
7250 ' begin: no ADC block 
7260 OUTPUT 707 USING 8800 ; V6,F6 

7270 3 1 =0 0 ON TIMERft 1,T7+!5M5 GOSUB 7550 ! read disable 
7200 ! end:no ADC block 
7290 T7=0 0 RETURN 

7 300 ' end: amplitude read setup procedure = = 
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7310 ! bsgin: ampl i tude read pr ocedur e-=======-=======-=== 

7320 ! in:Si,V0 out: v2,a2 
733 0 ! hold for read enable 
7340 IF SI THEN 7350 ELSE 7330 
7350 ! read enable 

73iii0 TRIGGER 723 ! measure xmtr voltage 
7370 ENTER 709 ) AS ! LIA magnitude output 
7390. 02=A5»V5 ! volts receiver output 
7390 ENTER 723 j VI ! xmtr voltage 
7400 A2=V2/V1 

7410 ! check for LIA overload 
7420 IF ASM THEN 7430 ELSE 7450 
7430 ' alarm -overload 

7440 PRINT "LIA OVERLOAD" S PRINT F6;A5;V1 S BEEP 55,2000 
7450 ! end: overload test 

74i0 DI3P USING 35G0 FS , ZO , 1 000^ V 1 , 1 000 » V2 
7470 IF S3 AND S9 AND I3<50 THEN 7480 ELSE 7530 
7490 ! graphics storage 
7490 13=13+1 

7500 IF S9^2 THEN X(I3)=Z0 ELSE X(I3)=F6 
7510 Y( I3)=1000^V2 ' mV 
7520 ! end graphics storage 
7530 RETURN 

7540 I end : amplitude read procedure = = = = = = = = '= = = = =- = = = = '- = " 
7550 ! b c- Q i n : r e a d enable p r c c c d u r e =■= = = =•■- = = = = = = = = = = = = = = •==■ -- 
7 560 I called by t i m e r 111 ■- r e s e t s read enable flag si 
7570 OFF TIMER# 1 
7590 SMM 
7590 RETURN 

7600 I end: read enable procedure======-===============-= 
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7610 ! begin: move manipulator down Z9 nm============= -========« 

7620 I input: :9 

7630 ENTER 720 ; NO ! counter ref 
7640 IF ABS(29X,5 THEN 7690 

7650 IF 29>0 THEN OUTPUT 709 ;"CLS1" S N3=l ELSE OUTPUT 709 ;"CLS 
0" 3 N3=-l 3 Z9=-(Z9/R8) 

7660 ON TIMERS 2 , I P ( Z 9/R2-T2 ) GOSUB 7700 ! set stop time 
7670 S1,S2=0 ! read off, inmotion 
7690 RETURN 

7690 ! end: move manipulator proced ure = = = = - = = " = = = = = = = = “ = 

7700 ! begin: manipulator stop pr ocedure = =^ = = = = = " = ” = = = ~~ = -' ■" = = -- = * 

7710 I input :0,29,n0,n3 

7720 ! output revised :0,step :8 

7730 IF N3>0 THEN OUTPUT 709 ;"0PN1" ELSE OUTPUT 709 ;"0PN0" 

7740 OFF TIMER# 2 

7750 ENTER 720 ; Ml 

7760 ! loop until ctr stops 

7770 ENTER 720 ; N2 

7790 IF N2=N1 THEN 7790 ELSE N1=N2 3 GOTO 7760 
7790 ! motion stop 
7300 Z9=N3*(N1-N0)/R0 
7810 Z0=Z0+ZS 

7920 DISP USING 3790 ; Z0,Z9 
7630 S2"0 1 flag motion stop 
7340 ; auto drive control 
7350 IF 36 THEN 7860 ELSE 7950 
7360 ! begin:adjust drive 
7870 ENTER 709 ; A5 
7330 IF A5>.95 THEN 7890 ELSE 7950 
7890 I LIA out too high 

7900 IF ASM. 12 THEN V 6 = FNR(V6m5) ELSE V 6 = FNR ( V 6 * . 7 ) 

7910 OUTPUT 707 USING 8810 ; V 6 
7920 PRINT "drive: ";V 6 ; "V" 

7930 WAIT I5*T3M3 
•’940 GOTO 7970 

7950 ON TIMER# 1,T7M5M5 GOSUB 7550 ! set delay to read e.iaole 
7960 77^0 
7'’70 RETURN 

7990 ‘ e n d : m a n i p Lt 1 a 1 0 r stop 3 ' o c e d u r e = ~ = = = ■ = ■ 

7990 ' beg:.n: manipulator relocate pr ocedure = = = - = = “ 

8000 CLEAR 

9010 IF 10=0 THEN DISP "po 3 ;t:on ";9 INPUT ZO 

8020 DISP "move to";0 INPUT Z4 

9030 Z9=Z4-Z0 ! mm 

3040 GOSUB 7610 I move :9 mm 

8050 S4=l 0 RETURN 

8060 ' end: manipulator relocate pro:edure- = = = = “ = = = "- " = = = "~ = -— = 
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8070 ! begin: title string justi f i er====-=====“======="“=--““ 

8080 I0=LEN(A$) 

8090 IF 10 MOD 32>0 THEN A$ = ATJ!" " 0 GOTO 8080 
8100 RETURN 

8110 i end title j ustifier " = = = = = = = = = = = " = === = = = "" = = = = = = = = = •-“ = •“= 

8120 ! ADC/:<ffltr toggle routine = -= = = = = = = = = = = = = = = = - = ”=- === = ---“- 

8130 IF 10=0 THEN 8140 ELSE 3170 

8140 ! operator call 

8150 S4 = l ! -fast re-fresh menu 

8160 CLEAR 0 DISP "1 = P}T. C^PIH 3-TlH 4 = cal";0 INPUT 13 
3170 ! process conunand 

8130 IF IB>0 THEN S6=l ELSE S6=0 0 I3=-I3 ! I8>0 turns ADC on 

8190 18=13 MOD 5 

3200 IF 13=0 THEN 3420 

8210 OUTPUT 707 ) "APOVO" 

3220 ON 13 SOTO 8230,3230,3330,3330 
8230 ' PIT 



3240 


S5 = 0 0 


ST = " 


P)T" 




8250 


OUTPUT 


706 


USING 


8820 


3260 


OUTPUT 


714 


USING 


3320 


8270 


GOTO 84 


20 






8230 


’ P)H 








S290 


55 = 0 0 


ST = " 


PIH" 




3300 


OUTPUT 


706 


USING 


3320 


5310 


OUTPUT 


714 


USING 


3820 


0320 


GOTO 34 


20 






3330 


' T3H 








3340 


35 = 1 0 


S$ = '' 


T)H" 




8350 


OUTPUT 


706 


USING 


3820 


3360 


OUTPUT 


714 


USING 


3320 


8370 


GOTO 84 


20 






3330 


! p r e a m 


p calibration 


5390 


35 = 0 0 




Plcal " 




3400 


OUTPUT 


706 


USING 


9820 


34 10 


OUTPUT 


714 


USING 


3820 


8420 


PRINT J 


SING 


“ ' H / " 


* 8-f 


5430 


13=0 0 


RETURN 




34 40 


1 end; 


:< fn t r 


/adc t 


0 g g 1 c 



'' B 1 A 1 " 
" 3 2 A 2 " 



"B2A1" 
" B 2 A 2 " 



"B2A2" ' 
”5303" 



"33Ar’ 
"BTAl " 
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8450 ! Images for formatted Qutput = = = = " = = = = = = = = = = -" = = = ""“ = " = "•“ = ■« 
8460 IMAGE /////K 

8470 IMAGE / / " t i me : " , 6D/ " cen t er freq measur ement " / " est i mat e: , 5D/ 
"range: " , 30/ " steps! : s : " , 30/ 



3430 


IMAGE 


#,"»",4D,2X,S.4DE 








8490 


IMAGE 


#,30.20 








3500 


IMAGE 


/"peak:", 50. 20, "He" 


/5X3D.5D, 


"mV"/5XS, 


3DE, "mV/Hz''2"/"xm 


tr : " , 


K, "mV" 










3510 


IMAGE 


/"resonance center 


frequency 


"/5D. 2D , " 


4. / - *« u~ " 


3520 


IMAGE 


/"1/2 power pt ob scur ed " / 40 , 


2(XK) 




8530 


IMAGE 


"low 1/2:", 50.0 








3540 


IMAGE 


"high 1/2:", 50.0 








3550 


IMAGE 


"resonance quality 


factor:". 


2D. 2D 




3560 


IMAGE 


/"search limit 3;<ceedcd",5D 






3570 


IMAGE 


/ " t : m 0 : " , 6 D / " Q in e a 5 u r 0 m e n t " / 


"center f 


req; " ,5D,2D/ '-depl 


h" ,40 


,0/ K 


/ 








8530 


IMAGE 


" 0 " ,4D,X4D.D,6D.2D, 


30.50 






3590 


IMAGE 


/K/"Fn-" ,K, "+/-" ,20 


. 20/"Cn=" 


? 1 ‘ 1 


2D. 22/ "Qn-'^ ,V/ 


Tl = ", 


i< , " d 3 " 


/"M(H)^=",K,"dB" 








3600 


IMAGE 


/"0-return to main 


l-LIA"/"2 


“water 


2 ■ - drive". " 


4 •" c 0 (n m ^ n t s 


3 -constants" 








2610 


IMAGE 


" 1 0 c k i n a n a 1 y c e r " / " 


input ",X 


2D. 4D, "mV 


" / " 0 p a n d " . X 2 C / *' t 



imo ",X50,"i<",20,"^",X50,"ms" 

3620 IMAGE "depth! m m ) : " , / " v o I u m e i m ■' 3 ) : " , 0 . 6 D 

2630 IMAGE K/"T 1 one t h = " , K , " mm " / " cap ac i t anc e ~ , K , " pF “ / " H lerq-!"-" 
,K,"mm"/ 

8640 IMAGE "water t emp er a t ur e " / 20 . 0 , " C" 

2630 IMAGE "6-titlo 7-plot c n/ of f " / " S^-p cs : t ; cn 

dcr specs" 

3660 IMA3E " , 40 . 20 , 2XK/ " ",2(XS.3DE) 



6670 


IMAGE 


p, " M ■» ” 


,K,7X,K, "#" ,20 












3690 


IMAGE 


3D. D/'m.n 


" ,K,6X ,K, " )t" , 2D 












36'^ 0 


IMAGE 


/:d," pt 


c"/"spd-" ,40. 20, " + 


.2D. 


*- " 1 


It 


m / z 


11 / •• _ 11 i. tl ^ <1 

' , . . Im . . 


3 70 0 


IMAGE 


K/"T cor 


rn: " ,320.20, "dB"/"H c 


0 r r n 


. M 

* 4 




n r\ 


2D. " d D " / 


3710 


IMAGE 


"pda", 3D 


.20,;(X.4DE) 












3720 


IMAGE 


/ " t i m 0 : " 


,60," m a ;; / m i n i o c a t o r 


11 / M ^ 
/ a 


t a r 


f 


. :i 

• 1 


4D , D , " r, t 00 .: ■ 


30.0/ 

2730 


IMAGE 


"mi .n / .ma:< 


: ,40. 0, ".nm"/3X3D. 50, 


"mV " 


./ II - 


1 


0 0 0 


; " . S . 4 D E / " - r 


. " c 

. , VJ . 


:de/ 
















3 7 4 0 


IMAGE 


" p r 0 ,;i m p 


gain: ",30.20, "dB § " 




mV" 








3750 


IMAGE 




,"V"/" It-" ,K, "mA"/" 


Eth- 


1 • • 


1 


"uV 


■1 


3760 


IMAGE 


/ / "-freq: 


" ,50. 20/ "sound spd:", 


4D. 2 


0/ " 


G! 


. >' 

• J 


3D. 3D 


3770 


IMAGE 


"volume: 


" ,K/ "density: ",40.20/ 


" Jn: 


",K 


/ 






8790 


IMAGE 


"T: " , 4D, 


2D,"rJB re lV/uPa"/"H: 


",4D 


. 20 


1 


"dB 


K .-1 1 W / 


3790 


IMAGE 


"•itop: " , 


40.20,50.40 












8900 


IMAGE 


" A P " 2D. 


50, "VG F.9",4D,"H2" 












8310 


IMAGE 


"AF " ’,2D. 


50, "VO" 
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9320 IMAGE #,4A 

9330 ! snd: i nages=============== 

8840 ! begin: end program routine 
3850 OUTPUT 709 ;"RS" 

8860 OUTPUT 707 ; "APOVO" 

8870 END 
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